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SUBSPACES 
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Abstract. A general classification of linear differential and finite-difference 
operators possessing a finite-dimensional invariant subspace with a polyno- 
mial basis (the generalized Bochner problem) is given. The main result is 
that any operator with the above property must have a representation as a 
polynomial element of the universal enveloping algebra of some algebra of 
differential (difference) operators in finite-dimensional representation plus an 
operator annihilating the finite-dimensional invariant subspace. In low dimen- 
sions a classification is given by algebras 3/2 (R) (for differential operators in 
R) and s?2(R)q (for finite-difference operators in R), osp(2, 2) (operators in 
one real and one Grassmann variable, or equivalcntly, 2x2 matrix operators in 
R), sZs(R), s/2(R)ffis^2(R) and gl'2(R) xR r+1 , r a natural number (operators 
in R 2 ). A classification of linear operators possessing infinitely many finite- 
dimensional invariant subspaces with a basis in polynomials is presented. A 
connection to the recently-discovered quasi-exactly-solvable spectral problems 
is discussed. 

S. Bochner (1929) asked about a classification of differential equations 

T<p - tip (0) 

where T is a linear differential operator of fcth order in one real variable x G R and 
e is the spectral parameter, having an infinite sequence of orthogonal polynomial 
solutions |l| (see also 0). 

Definition 0.1. Let us give the name of the generalized Bochner problem to the 
problem of classification of linear differential (difference) operators, for which the 
eigenvalue problem (0) has a certain number of eigenfunctions in the form of a 
finite-order polynomial in some variables. 

Following this definition the original Bochner problem is simply a particular case. 
In || a general method has been formulated for generating eigenvalue problems for 
linear differential operators, linear matrix differential operators and linear finite- 
difference operators in one and several variables possessing polynomial solutions. 
The method was based on considering the eigenvalue problem for the representation 
of a polynomial element of the universal enveloping algebra of the Lie algebra in a 
finite-dimensional, 'projectivized' representation of this Lie algebra q|. Below it is 
shown that this method provides both necessary and sufficient conditions for the 
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existence of polynomial solutions of linear differential equations and a certain class 
of finite-difference equations. 

The generalized Bochner problem can be subdivided into two parts: (i) a clas- 
sification of linear operators possessing an invariant subspace (subspaces) with a 
basis in polynomials of finite degree and (ii) a description of the conditions under 
which such operators are symmetrical. This paper will be devoted to a solution 
of the first problem; as for the second one, only the first step has been done (see 
below) . The plan of the paper is the following: Section 1 is devoted to the case of 
differential operators in one real variable; in Section 2 finite-difference operators in 
R are treated; operators in one real and one Grassmann variables are considered 
in Section 3, while the case of operators in R 2 is given in Section 4. The general 
situation is described in the Conclusion. 

1. Ordinary differential equations 
Consider the space of all polynomials of order n 

V n = (l,X,X 2 ,...,X n ), (f) 

where n is a non-negative integer and x G R. 

Definition 1.1. Let us name a linear differential operator of the A:th order, 
quasi- exactly- solvable, if it preserves the space V n . Correspondingly, the operator 
Ek , which preserves the infinite flag Vo C V\ C V% C • • • C V n C . . . of spaces of 
all polynomials, is named exactly-solvable. 

Lemma 1.1. (i) Suppose n > (fc — 1). Any quasi- exactly- solvable operator Tk, can 
be represented by a k-th degree polynomial of the operators 

J + = x 2 d x — nx, 
J° = xd x -2, (2) 
J~ = d x , 

(the operators (2) obey the s^(R) commutation relations [] ). If n < (k — I), the 
part of the quasi-exactly-solvable operator 7\. containing derivatives up to order n 
can be represented by an nth degree polynomial in the generators (2). 

(ii) Conversely, any polynomial in (2) is quasi-exactly solvable. 

(Hi) Among quasi-exactly-solvable operators there exist exactly- solvable operators 
EkCT k . 

Comment 1. If we define the universal enveloping algebra U g of a Lie algebra 
g as the algebra of all ordered polynomials in generators, then T& at k < n + 1 is 
simply an element of the universal enveloping algebra C/ s ; 2 (r) of the algebra s?2(R) 
taken in representation (2). If k > n + 1, then Tk is represented as an element of 
U s i 2 (r) plus B £ X n+i i where B is any linear differential operator of order not higher 
than (k — n — 1). 

Proof. The essential part of the proof is based on the Burnside theorem (see, e.g., 

m 

x The representation (2) is one of the 'projectivized' representations (see ||). This realization 
of si2(R-) has been derived at the first time by Sophus Lie. 

2 Burnside theorem is a particular case of more general Jacobson theorem (see JiJ, Chapter 
XVII. 3). I am grateful to V. Kac for this comment. 
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Let A\ , . . . , Ak be linear operators in a linear space E over real numbers, 
dimE < oo. Let us assume that there is no linear space L over real 
numbers, < dimL < dimE , such that A; : L i— > L for all i 
1.2.... , fc. Then any linear operator acting in E can be represented as 
a polynomial in Ai , . . . , A^. 
The operators J '° act in P n irreducibly and therefore the theorem can be applied. 
Also there exist operators B having V n as a kernel, B : V n >— > 0. Clearly, those 
operators have a form B(x, d x ) ^ xn+1 , where B(x,d x ) is any linear differential op- 
erator, and they make no contribution in Tk if k < n+ 1. It completes the proof of 
part (i). Parts (ii) and (iii) are easy to prove based on part (i).[] □ 

Since s^(R-) is a graded algebra, let us introduce the grading of generators (2): 
deg(J+) = +1 , deg{J°) = , deg(J~) = -1, (3) 

hence 

degl(J + r+(J°r o (J-T-} = n + -n_. (4) 
The grading allows us to classify the operators in the Lie-algebraic sense. 

Lemma 1.2. A quasi- exactly- solvable operator C U s i 2 ^ has no terms of posi- 
tive grading, if and only if it is an exactly-solvable operator. 



Theorem 1.1. Let n be a non-negative integer. Take the eigenvalue problem for a 
linear differential operator of the kth order in one variable 

T k ip = Eip , (5) 

where Tfc is symmetric. The problem (5) has (n + 1) linearly independent eigen- 
functions in the form of a polynomial in variable x of order not higher than n, if 
and only if Tk is quasi-exactly-solvable. The problem (5) has an infinite sequence 
of polynomial eigenf unctions, if and only if the operator is exactly-solvable. 

Comment 2. The " if " part of the first statement is obvious. The " only if " part 
is a direct corollary of Lemma 1.1 . 

This theorem gives a general classification of differential equations 

k 

^2aj{x)(p {l) (x) = eip{x) (6) 
j=o 

having at least one polynomial solution in x. The coefficient functions aj(x) must 
have the form 

k+j 

i=0 

The explicit expressions (7) for coefficient function in (6) are obtained by the substi- 
tution (2) into a general, fcth degree polynomial clement of the universal enveloping 
algebra tZ^^R). Thus the coefficients Oj^ can be expressed through the coefficients 
of the fcth degree polynomial clement of the universal enveloping algebra U s i 2 ^y 
The number of free parameters of the polynomial solutions is defined by the num- 
ber of parameters characterizing a general fcth degree polynomial element of the 



5 A detailed, technical proof is presented in Jl3| . 



4 



ALEXANDER TURBINER 



universal enveloping algebra C/ s ; 2 (r). A rather straightforward calculation leads to 
the following formula 

par{T k ) = {k + lf (8) 

where we denote the number of free parameters of operator Tk by the symbol 
par(Tk). For the case of an infinite sequence of polynomial solutions expression (7) 
simplifies to 

j 

in agreement with the results of H.L. Krall's classification theorem ^)(see also 0). 
In this case the number of free parameters is equal to 

, , (fc + l)(A; + 2) , s 

par(E k ) = ± >j i (10) 

In the present approach Krall's theorem is simply a description of differential oper- 
ators of /cth order in one variable preserving a finite flag Vo C V\ C V% C ■ ■ • C Vu 
of spaces of polynomials. One can easily show that the preservation of such a set 
of polynomial spaces implies the preservation of an infinite flag of such spaces. 

One may ask a more general question: which non-degenerate linear differential 
operators have a finite-dimensional invariant sub-space of the form 

(a(z),a(z)(3(z),...,a(z)[3(zr) , (11) 

where a(z),(3(z) are some functions. Such operators are easily obtained from the 
quasi-exactly-solvable operators (see Lemma 1.1) making the change of variable 
x — (3{z) and the "gauge" transformation T = a(z)Ta(z) _1 . Since any one- or 
two-dimensional invariant sub-space can be presented in the form (11), the impor- 
tant general statement takes place: there are no linear operators possessing one- 
or two-dimensional invariant sub- space with an explicit basis other than given by 
Lemma 1.1 . Therefore for any linear operator possessing one- or two-dimensional 
invariant sub-space with an explicit basis the eigenvalue problem (5) can be reduced 
to the form (16) or (15), respectively (see below). 

Let us consider the second-order differential equation (6), which can possess poly- 
nomial solutions. From Theorem 1.1 it follows that the corresponding differential 
operator must be quasi-exactly-solvable and can be represented as 

T 2 = C++J+J+ + c +0 J + J° + C+-J+J- + c -J°J- + C—J-J-+ 

c+J+ +cqJ° + c-J- +c, (12) 

where c a p,c a ,c € R. The number of free parameters is pariT^) = 9. Under the 
condition c ++ = c + o = c + = 0, the operator T2 becomes exactly-solvable (see 
Lemma 1.2) and the number of free parameters is reduced to par{E2) = 6. 

Lemma 1.3. If the operator (12) is such that 

n 

c ++ = and c+ = (— — m)c+o 7 a t some m = 0, 1, 2, . . . (13) 

then the operator T2 preserves both V n cmd V m - In this case the number of free 
parameters is par(T 2 ) = 7. 
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In fact, Lemma 1.3 means that T 2 ( J Q (n), c a p, c a ) can be rewritten as T 2 (J a (m), 
c q/3' c q)- A s a consequence of Lemma 1.3 and Theorem 1.1, in general, among 
polynomial solutions of (6) there are polynomials of order n and order m. 

Remark. From the Lie-algebraic point of view Lemma 1.3 means the existence 
of representations of second-degree polynomials in the generators (2) possessing two 
invariant sub-spaces. In general, if n in (2) is a non-negative integer, then among 
representations of fcth degree polynomials in the generators (2), lying in the uni- 
versal enveloping algebra, there exist representations possessing 1,2, k invariant 
sub-spaces. Even starting from an infinite-dimensional representation of the origi- 
nal algebra (n in (2) is not a non-negative integer), one can construct the elements 
of the universal enveloping algebra having finite-dimensional representation (e.g., 
the parameter n in (13) is non-integer, however T 2 has the invariant sub-space of 
dimension (m + 1)). Also this property implies the existence of representations of 
the polynomial elements of the universal enveloping algebra, which can be obtained 
starting from different representations of the original algebra. 

Substituting (2) into (12) and then into (6), we obtain 

- P 4 {x)d 2 Mx) + Pz{x)d x y{x) + P 2 (x)<p(x) = e<p(x), (14) 

where the Pj(x) are polynomials of jth order with coefficients related to c a f3,c a 
and n. In general, problem (14) has {n + 1) polynomial solutions. If n = 1, as a 
consequence of Lemma 1.1, a more general spectral problem than (14) arises 

- F 3 (x)d%<p(x) + Q 2 {x)d x ip{x) + Qi{x)ip(x) = £<p(x), (15) 

where F$ is an arbitrary real function of x and Qj{x),j = 1,2 are polynomials of 
order j, possessing only two polynomial solutions of the form (ax + b). For the case 
n = (one polynomial solution) the spectral problem (6) becomes 

-F 2 {x)dly{x) + F 1 (x)d s <p(x) + Qo<p(x) = E (p(x), (16) 

where ^2,1(2;) are arbitrary real functions and Qo is a real constant. After the 
transformation 

tt(z) = *{x{z))e- A (*\ (17) 

where x — x{z) is a change of the variable and A(z) is a certain real function, one 
can reduce (14)-(16) to the Sturm-Liouville-type problem 

{-d 2 z + V(z))9(z) = e*(«), (18) 

with the potential, which is equal to 

V(z) = (A') 2 - A" + P 2 (x(z)) , 

if 




for the case of (14). If the functions (17), obtained after transformation, belong 
to the £2(f)-spaceP| , we reproduce the recently discovered quasi-exactly-solvable 
problems @], where a finite number of eigenstates was found algebraically. For 
example, 

T 2 = -4J° J- + 4aJ+ + 4bJ° - 2(n + 1 + 2k) J~ + b(2n + 1 + 2fc) (19) 



4 Depending on the change of variable x = x(z), the space T> can be whole real line, semi-line 
and a finite interval. 
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leads to the spectral problem (18) (x = z 2 ) with the potential 

V{z) = a 2 z 6 + 2abz 4 + [b 2 - (4n + 3 + 2k)a]z 2 , (20) 

for which at k = 0(k = 1) the first (n + 1) eigenfunctions, even (odd) in x, can be 
found algebraically. Of course, the number of those 'algebraized' eigenfunctions is 
nothing but the dimension of the irreducible representation (1). 

Taking different exactly-solvable operators E2 for the eigenvalue problem (6) one 
can reproduce the equations having the Hermite, Laguerre, Legendre and Jacobi 
polynomials as solutions It is worth noting that for general exactly-solvable 
operator Ei the eigenvalues are quadratic in number of eigenstate and can be 
presented as follows 

e« = coon 2 + con + const 

(for details see ||). 

Under special choices of the general element Ei(Ee,Eg), one can reproduce all 
known fourth- (sixth-, eighth-)order differential equations giving rise to infinite se- 
quences of orthogonal polynomials (see Q and other papers in this volume) . 

Recently, A. Gonzalez-Lopez, N. Kamran and P. Olver |)| gave the complete 
description of second-order polynomial elements of U s i 2 (j^ leading to the square- 
integrablc eigenfunctions of the Sturm-Liouville problem (18) after transformation 
(17). Consequently, for second-order ordinary differential equations (14) the com- 
bination of their result and Theorem 1.1 gives a general solution of the Bochncr 
problem as well as the more general problem of classification of equations possessing 
a finite number of orthogonal polynomial solutions. 

2. Finite-difference equations in one variable 

Let us introduce the finite-difference analogue of the differential operators (2) 


J+ = x 2 D - {n}x 

J° = xD - h (21) 
J- = D, 

where h = ^ {271+2} ^ > i. n } = Y-q ^ s * ne q uan t um symbol (or g-number), q G R is 
a number characterizing the deformation, Df(x) — ^jlff^ " is a shift or a finite- 
difference operator (or the so-called Jackson symbol (sec flljl)). The operators (21) 
after multiplication by some factors 

~ _ g-" {2n + 2} J0 
q+1 {n+l} 
j± = q -n/2J± 

( see [fio|| ) form a quantum si2(R) 9 algebra with the following commutation relations 

qfr - r f = -r 

q 2 ] + r - r~j + = -(<Z+l)j° (22) 



5 For instance, setting the parameter a = in (19), the equation (15) converts to the Hermite 
equation (after some substitution). 
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(this algebra corresponds to the second Witten quantum deformation of s?2 in the 
classification of C. Zachos |l2|]). If q — > 1, the commutation relations (22) reduce 
to the standard s^R) ones. A remarkable property of generators (21) is that, if 
n is a non-negative integer, they form the finite-dimensional representation. For q 
others than root of unity this representation is irreducible. 

Similarly as for differential operators one can introduce quasi-exactly-solvable 
Tfc and exactly-solvable operators Ek- 

Lemma 2.1. (i) Suppose n > [k— 1). Any quasi- exactly- solvable operator Tk, can 
be represented by a kth degree polynomial of the operators (21). If n < (k — 1), the 
part of the quasi-exactly-solvable operator Tk containing derivatives up to order n 
can be represented by a nth degree polynomial in the generators (21). 

(ii) Conversely, any polynomial in (21) is quasi-exactly solvable. 

(Hi) Among quasi-exactly-solvable operators there exist exactly-solvable operators 
E k df k . 

Comment 3. If we define an analogue of the universal enveloping algebra U g for the 
quantum algebra g as an algebra of all ordered polynomials in generators, then a 
quasi-exactly-solvable operator Tk at k < n+ 1 is simply an element of the 'universal 
enveloping algebra' f/ s / 2 (R,) of the algebra s^R)? taken in representation (21). If 
k > n + 1, then Tk is represented as an element of U s i 2 ^ q plus BD n+1 , where B 
is any linear difference operator of order not higher than (k — n — I). 

Similar to s^R) , one can introduce the grading of generators (21) of sl2(R) g 
(see (3)) and, hence, of monomials of the universal enveloping J7 s z 2 (r), i (see (4)). 

Lemma 2.2. A quasi-exactly-solvable operator Tk C U s i 2 ^ has no terms of pos- 
itive grading, iff it is an exactly-solvable operator. 

Theorem 2.1. Let n be a non-negative integer. Take the eigenvalue problem for a 
linear difference operator of the k-th order in one variable 

T k ip(x) = eip(x), (23) 

where T^ is symmetric. The problem (23) has (n + 1) linearly independent eigen- 
functions in the form of a polynomial in variable x of order not higher than n, if 
and only ifT^ is quasi-exactly-solvable. The problem (23) has an infinite sequence 
of polynomial eigenf unctions, if and only if the operator is exactly- solvable E^. 

Comment 4- Saying the operator Tj. is symmetric, we imply that, considering 
the action of this operator on a space of polynomials of degree not higher than 
n, one can introduce a positively-defined scalar product, and the operator Tk is 
symmetric with respect to it. 

This theorem gives a general classification of finite-difference equations 

k 

J2 aj {x)D j ip{x) = e<p(x) (24) 
j=o 

having polynomial solutions in x. The coefficient functions must have the form 

k+j 

a j( x ) = ^Z^j,^ 1 - ( 25 ) 

i=0 
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In particular, this form occurs after substitution (21) into a general fcth degree 
polynomial element of the universal enveloping algebra U s [ 2 (rn . It guarantees the 
existence of at least a finite number of polynomial solutions. The coefficients ctjj 
are related to the coefficients of the fcth degree polynomial element of the univer- 
sal enveloping algebra ?7 s / 2 (r,) . The number of free parameters of the polynomial 
solutions is defined by the number of free parameters of a general k-th order polyno- 
mial element of the universal enveloping algebra U s i 2 ^ q A rather straightforward 
calculation leads to the following formula 

par(f k ) = (k + l) 2 + 1 

(for the second-order finite-difference equation par(T 2 ) = 10). For the case of an 
infinite sequence of polynomial solutions the formula (25) simplifies to 



= ^2°>3,i xi ( 26 ) 



=0 



and the number of free parameters is given by 

, - . (k + l)(k + 2) 
par{E k ) = K - >j + 1 

(for k = 2, par(E 2 ) — 7). The increase in the number of free parameters com- 
pared to ordinary differential equations is due to the presence of the deformation 
parameter q. In || one can find a description in the present approach of the q- 
deformed Hermite, Laguerre, Legendre and Jacobi polynomials (for definitions of 



these polynomials see [11 



Lemma 2.3. If the operator (see (12)) is such that 

c ++ = and c + = (n — {m})c+o , at some m = 0, 1, 2, . . . (27) 

then the operator T 2 preserves both V n and V m , and polynomial solutions in x with 
8 free parameters occur. 

As usual in quantum algebras, a rather outstanding situation occurs if the 
deformation parameter q is equal to the root of unity. For instance, the follow- 
ing statement holds. 

Lemma 2.4. If a quasi- exactly- solvable operator T k preserves the space V n and the 
parameter q satisfies to the equation 

q n = 1 , (28) 

then the operator Tk preserves an infinite flag of polynomial spaces Vq C V n C 
P21.C-C Vkn C .... 

It is worth emphasizing that, in the limit as q tends to one, Lemmas 2.1,2,3 and 
Theorem 2.1 coincide with Lemmas 1.1,2,3 and Theorem 1.1, respectively. Thus 
the case of differential equations in one variable can be treated as particular case 



6 For quantum sl2(R) q algebra there are no polynomial Casimir operators (see, e.g., Jl^]). 
However, in the representation (21) the relationship between generators analogous to the quadratic 
Casimir operator 

qj+j- - J°J° + ({n + 1} - 2h)J° = h{h -{n + 1}) 
appears. It reduces the number of independent parameters of the second-order polynomial element 
of f/ s i 2 (R) (! . It becomes the standard Casimir operator at q — » 1. 
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of finite-difference ones. Evidently, one can consider the finite-difference operators, 
which are a mixture of generators (21) with the same value of n and different g's. 

3. Operators in one real and one Grassmann variable 

Define the following space of polynomials in x, 9 

V NM = (x°,x 1 ,...,x N ,x°9,x 1 9,...,x M 9) (29) 

where N, M are non-negative integers, x G R and 9 is a Grassmann (anticommut- 
ing) variable. 

The projectivized representation of the algebra osp(2, 2) is given as follows. 
The algebra osp(2, 2) is characterized by four bosonic generators T ±,0 5 J and 
four fermionic generators Qi,2, Q\ 2 and given by the commutation relations 

[T°,T ± ]=±T ± , [T + ,T~] = -2T , [J,T Q ]=0 ,a = +,-,0 
{Qi-,Q 2 } = -T- , {Q 2! Qi}-T+ , 
{Q 2 ,Q 2 }) = +J 

{Qi,Qi} = {Q 2 ,Q 2 } = {Qi,Q 2 } = , 

(Qi,Qi} = {Q 2 ,Q 2 } = {Qi,Q 2 } = , 
, [Q 2 ,T+] = , [Q u T-]=0 , [Q 2 ,T-} = -Q! 



\{{Qx,Qi] 



1({Qi,Qi}-{Q 2 ,Q2}) = t° 



[Qi,T + ] 



Q2 





[g 2 ,r + ] = -Qi , [Qi,t-}=' 



[Q 2 ,T-} = 



[Qi, 2) T°} = ±Iq 1:2 



\Qi,2,T )=tIQi, 2 



[Qi,a,J] 



1 



h,2 , [Qi, 2 ,J] 



1- 



;Qi, 



(30) 



This algebra has the algebra sl 2 (R) ® R as sub-algebra. 

The algebra (30) possesses the projectivized representation 

T+ = x 2 d x -nx + x9d e , 



rpO 



(31) 



9.r 



2 2 



for bosonic (even) generators and 

Q = 



~Q{ 




~dg- 


, Q = 


Q{ 




x9d x — n9 


_Q 2 _ 




xdg_ 


_Q 2 _ 




-9d x 



(32) 



for fermionic (odd) generators, where a: is a real variable and 9 is a Grassmann 
one. Inspection of the generators shows that if n is a non-negative integer, the 
representation (31), (32) is a finite-dimensional representation of dimension (2n+l). 
The polynomial space V n .n-i describes the corresponding invariant sub-space. 

Definition 3.1. Let us name a linear differential operator of the fc-th order, Tk(x, 9) 
quasi-exactly-solvable if it preserves the space V n , n -i- Correspondingly, the 
operator Ek(x, 9) € Tk(x, 9), which preserves the infinite flag 'Po.o C 7^1,0 C V 2 ,i C 
■ • ■ C Vn.n-i C . . . of spaces of all polynomials, is named exactly-solvable. 
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Lemma 3.1. Consider the space V n , n -i- 

(i) Suppose n > (k — 1). Any quasi-exactly-solvable operator Tk(x,0), can be 
represented by a kth degree polynomial of the operators (31), (32). If n < (k — 1), 
the part of the quasi-exactly-solvable operator Tk(x,6) containing derivatives in x 
up to order n can be represented by an nth degree polynomial in the generators (31 ), 
(32). 

(ii) Conversely, any polynomial in (31), (32) is a quasi-exactly solvable operator, 
(in) Among quasi-exactly-solvable operators there exist exactly-solvable operators 

E fe C T k (x,0). 

Let us introduce the grading of the bosonic generators (31) 

deg{T+) = +1 , deg{J, T°) = , deg(J~) = -1 (33) 
and fermionic generators (32) 

deg(Q 2 ,Q 1 ) = +^ , deg(Q u Q 2 ) = ~ (34) 
Hence the grading of monomials of the generators (31), (32) is equal to 
deg[(T+) n + (T°) no (J) w (T-) n - Q 1 mi Q 2 m2 Q 1 Wl Q 2 ^ 2 } = 

(n+ — n_) — (mi — m 2 — mi + m2)/2 (35) 

The n's can be arbitrary non-negative integers, while the m's are equal to cither 
or 1. The notion of grading allows us to classify the operators T k (x,9) in the 
Lie-algebraic sense. 

Lemma 3.2. A quasi-exactly-solvable operator T k (x, 9) C U osp ^ 2 ^ has no terms 
of positive grading other than monomials of grading +1/2 containing the generator 
Qi or Q 2 , iff it is an exactly- solvable operator. 

Theorem 3.1. Let n be a non-negative integer. Take the eigenvalue problem for a 
linear differential operator in one real and one Grassmann variable 

T k {x,6)y = e<p (36) 

where T& is symmetric. In general, the problem (36) has (2n + 1) linearly indepen- 
dent eigenf unctions in the form of polynomials in variables x, 9 of order not higher 
than n, if and only ifT k is quasi-exactly-solvable. The problem (36) has an infinite 
sequence of polynomial eigenf unctions, if and only if the operator is exactly-solvable. 

This theorem gives a general classification of differential equations 

i=k,j=l 

J2 aij{x,0)<p*f{x,6) = e<p(x,e) , (37) 

i,j=0 

where the notation <p x 'g means the ith order derivative with respect to x and jth 
order derivative with respect to 9, having at least one polynomial solution in x, 9, 
thus resolving the generalized Bochner problem. Suppose that k > 0, then the 
coefficient functions aij(x,9) should have the form 

k+i k+i— 1 

a-i,o{x,0) = ^a ij0 , p x p + 9 ^ a ifi ^ x p 

p=0 p=0 
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k+i-1 k+i-1 

a it i(x,6) = ^ a iA,P xP + 9 xP ( 38 ) 

p=0 p=0 

The explicit expressions (38) are obtained by substituting (31), (32) into a general, 
the fcth order, polynomial element of the universal enveloping algebra U osp ( 2 , 2 y 
Thus the coefficients aij jP can be expressed through the coefficients of the fc-th 
order polynomial clement of universal enveloping algebra f/ «p(2,2)- The number of 
free parameters of the polynomial solutions is defined by the number of parameters 
characterizing a general, kth order polynomial element of the universal enveloping 
algebra U osp ( 2 .2)- However, in counting parameters some relations between genera- 
tors should be taken into account, specific for the given representation (31), (32), 
like 

2T+J - Q X Q 2 = nT+ , 
T+Q 1 - T°Q 2 = -Q 2 , 

T+Q 2 + T'Q, = Q t , 

71 

JQi = 2 ^1 ' 

T+T- - T°T° - J J + T° = -2(„_i) , 

J J = (n + i)J - J(n + 1) , 

Q1Q1 + Q2Q2 - 2n J = -n(n + 1) , 
2T°J + Q 1 Q 1 - (n + l)T° - nJ = -|(n+l), 

T-Q 2 - T°Q 1 - (| + l)Qi , 
T~Q X - T°Q 2 = ^y^)Q 2 , 

JQi = \ Qi , 

2JT- - = (n + l)T~ (39) 

between quadratic expressions in generators (and the ideals generated by them). 
Straightforward analysis leads to the following formula for the number of free 
parameters 

par(T k (x,0)) = 4fc(fc + 1) + 1 , (40) 

For the case of an exactly-solvable operator (an infinite sequence of polynomial 
solutions of Eq. (37)), the expressions (38) simplify and reduce to 

i i—1 

a it0 (x,8) = ^a ifi ^ p x p + 0^a lfi , p x p 

p=0 p=0 
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i i—X 

a it i(x,6) = aj,i, p x p + ai,i,p% p (41) 

p=0 p=0 

Correspondingly, the number of free parameters reduces to 

par(E fe (x,0)) = 2fc(fc + 2) + 1 (42) 

Hence, Eq. (37) with the coefficient functions (41) gives a general form of eigenvalue 
problem for the operator which can lead to an infinite family of orthogonal 
polynomials as eigenfunctions. If in (41) we put formally all coefficients, a^o.p and 
ai t \{x,ff) equal to zero, we reproduce the eigenvalue problem for the differential 
operators in one real variable, which gives rise to all known families of orthogonal 
polynomials in one real variable (see (l3|), 

3.1. Second-order differential equations in x,9. Now let us consider in more 
detail the second-order differential equation Eq. (37) , which can possess polynomial 
solutions. As follows from Theorem 3.1, the corresponding differential operator 
T 2 (x, 9) should be quasi-exactly-solvable. Hence, this operator can be expressed in 
terms of osp(2, 2) generators taking into account the relations (39) 

T 2 = C++T+T+ + c +0 T + r° + c+.T+T" + c _T°T- + c__T-T~+ 
c+jT+J + CojT°J + c-jT-J+ 
c +T T+Q 1 + c +2 T+Q 2 + c+iT+Qi + c +1 T+Q 2 + coiT°Qi+ 
c 2 T °Q2 + c-iT'Qx + c_zT-Q 2 + 
c+T+ + c T + c_T- + cjj + dQ 1 + c 2 Q 2 + c T Q 1 + Cl Q 2 + c (43) 
where c a p, c ai c are parameters. Following Lemma 3.2, under the conditions 

C++ = C+o = C +T = C +I = Cy = C+2 = C+j = C+ = , (44) 

the operator T2(x,9) becomes exactly-solvable (see (41)). 

Now we proceed to the detailed analysis of the quasi-exactly-solvable operator 
T 2 (x,9). Set 

c++ = (45) 

in Eq. (43). The remainder possesses an exceptionally rich structure. The whole 
situation can be subdivided into three cases 

c+2 ^ , c +T = (case I) (46) 
c+2 = , c +T ^ (case II) (47) 

c+2 = , c +T = (case III) (48) 

We emphasize that we keep the parameter n in the representation (31), (32) as a 
fixed, non-negative integer. 
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Case I. The conditions (45) and (46) are fulfilled (see Fig. 3.1). 
Case 1. 1.1. If 

(n + 2)c +0 + nc +J + 2c + = , 

c +2 = C T = ' 

(n+l)c 2 + 2c2 = , (49) 

then T 2 (x, 0) preserves V n , n -i and V n +i,n-i- 
Case 1.1.2. If 

(n + 4 + 2ra)c+o + nc + j + 2c + = , 

c +2 = C T = ^ ' 
C 02 = C2 = C -2 = . ( 50 ) 

at a certain integer m > , then T 2 (x, 0) preserves P n ,n-i and 7 , n +2+m,n-i- If to 
is non-integer, then T 2 {x, 9) preserves P n ,n-i and V<x>,n-i- 
Case 1.1.3. If 

(n + l)c +J + 2c + = , 
c+o = , 

C +2 = °T = i 
C 02 = C 2 = C -2 = . ( 51 ) 

then T 2 (x, 9) preserves the infinite flag of polynomial spaces the ? n + m ,n-i, 
to = 0,1, 2,.... 
Case J. 2. J. If 

(n - 3)c +0 + (n + l)c+j + 2c+ = , 
(n- l)c +5 = c T , 

(n-l)c 2 + 2^=0 , (52) 

then T 2 (x,9) preserves P^n-i and V n . n - 2 . 
Case 1.2.2. If 

3c +0 - c +J = , 
(2k + 2n + 4)c +0 + 2c+ = , 

C +2 = °T = ' 

(2fc - n + 3)c 02 - + 2c2 = , (53) 

at a certain integer k > 0, then T 2 (a;,0) preserves V n , n -i and ~Pk+2,k- 
Case 1.2.3. If 

c+o = c+,7 = c + = , 

c +2 = C T = ! 

C 02 = °2 = ' ( 54 ) 

then T 2 (x,9) preserves 'Pn.n— l and the infinite flag of the polynomial spaces T ? k+2.k? 
k = 0, 1, 2, ... . Note in general for this case c_2 ^= 0. 
Case J.3.J. If 

(n - 5 - 2m)c +0 + (n + l)c +J + 2c+ = , 

c +2 = C T = ! 
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c 02 — c 2 ~ c -2 — , (55) 

at a certain integer < m < (n — 3) , then T 2 (x, 9) preserves V n , n -i and V n , n -3-m- 
Case 1.3.2. If 

c+o = , 
(n + l)c +J + 2c+ = , 
c + 2 = C T = ! 

C 02 = C 2 = C -2 = > ( 56 ) 

then T2 (a;, 0) preserves V n ,n-\ and the sequence of the polynomial spaces V n ,n-3-m, 
to = 0,1,2,... ,(n- 3). 
Case 7.5.5. If 

(2fc + 1 - n)c +0 + (n + l)c+j + 2c+ = , 
(2m + 5)c +0 - c + ,7 = , 

C +2 = ^ = 1 
C 02 = C 2 = C -2 = . ( 57 ) 

at certain integers k,m > , then T 2 (x, #) preserves Pn.n-i and Vk+z+m.k- 
Case I.3.4. If 

c+o = c+J = c+ = , 

C +2 = ^ = ' 



C 02 ~ C 2 ~ C -2 



, (58) 

then T 2 (x, 9) preserves V n , n -i and the infinite flag of polynomial spaces Vk+3+ m ,k, 
k,m = 0, 1,2, ... (cf. Cases 1.1.3 and 1.2.3). 

Case II. The conditions (45) and (47) are fulfilled (see Fig. 3.II). 

Case II. 1.1. If 

(n + l)c +0 + (n + l)c +J + 2c+ = , 

c 2 = 0, (59) 

then T 2 (x,6) preserves V n ,n-i and 
Case II. 1.2. If 

(n + 3)c +0 + (n + l)c +J + 2c+ = , 
(n + 2)coi + 2ci = , 

c+i = c 2 = , (60) 

then T 2 (x, 6*) preserves V n , n -i and P ni „+i. 
Case II. 1.3. If 

(2fe + 5 + n)c +0 + (n + l)c +J + 2c+ = , 
coi = ci = , 
c+i = c 2 = , 

c_i=0, (61) 
at a certain integer k > , then T 2 (x, 9) preserves V n ,n-i and 7 : 'n,n+2+fc- 
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Case II.I.4. If 

c+o = 0, 
(n+l)c +J + 2c + = , 
coi = ci = , 
c+i = c 2 = , 

c-i = 0, (62) 

then T 2 (a:,#) preserves V n , n -i and the infinite flag of polynomial spaces V n , n +k, 
k = 0, 1, 2, . . . (cf. Cases 7.7.3, J. 2. 3 and 7.3.^). 
Case 77.2.7. If 

(n - 2)c +0 + nc+j + 2c + = , 

c+i = c 2 , (63) 

then T 2 (x,9) preserves V n ,n-i and V n -i,n-i- 
Case II.2.2. If 

(n + l)c +0 - c+ = , 
3c +0 + c+j = , 
nc i + 2ci = , 

c+i = c 2 = , (64) 

then T 2 (x,6) preserves V n , n -i and P n -i,n- 
Case II. 2.3. If 

(n - 2)c +0 + nc + j + 2c + = , 
(2k + l)c +0 + c+j = , 
c+i = c 2 = , 

coi = ci = c_i = , (65) 

at a certain integer k > 0, then T 2 (a:,#) preserves Pn^-i and P n -i.ii+t+i- 
Case II.24. If 

c+o = c+j = c + = , 
c+i = c 2 = , 

coi = ci = c_i = , (66) 

then T 2 (x, 6 1 ) preserves V n , n -i and the infinite flag of the polynomial spaces V n -i. n +k, 
k = 0, 1, 2, . . . (cf. Cases 1.7.3, 12.3, 7.3.^ and 77. 7. 4). 
Case 77.3.7. If 

(n - 4)c +0 + nc + j + 2c + = , 
(n - 2)cqi + 2 Cl = , 

c+i = c 2 = , (67) 

then T 2 (x,9) preserves V n ,n-i and V„-2,n-i- 
Case II.4.I. If 

(m — 2n)c + o + c + = , 
3c +0 + c + j = , 
(2m + 2 - n)c i + 2ci = , 

c+i = c 2 = , (68) 
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at a certain integer to > 0, then T 2 (x,9) preserves V n , n -i and V m ,m+\- 
Case 114.2. If 

(2m — n)c +0 + nc +J + 2c + = , 
(2k + 5)c +0 + 2c +J = , 
c+i = c 2 = 0, 

coi = ci = c_i = (69) 

at certain integers fc > 0, m > 0, then T 2 (x,#) preserves V n ,n-i and V m ,m+2+k- 
Case If 

c+o = c+j = c + = , 
c+i = c 2 = 0, 

coi = ci = c_i = (70) 

then T 2 (x,0) preserves l^n,n—i &nd the infinite flag of the polynomial spaces 
"Pm,m+i+fe, m, fe = 0,1,2,... (cf. Cases 1.1.3, 1.2.3, 1.3.4, H-l-4 and 11.24). 

Case III. The conditions (45) and (48) arc fulfilled (see Fig. 3.III). 

Case III. 1.1. If 

(2m - n)c +0 + nc+j + 2c + = , 
c+o + c+j = , 

(m - n)c+i + c 2 = 0, (71) 
at a certain integer to > 0, then T 2 (:r,0) preserves 1~*n,n—i &nd 7^m,m- 
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Case III. 1.2. If 

c+o = c+,7 = c+ = , 

c+i = c 2 = 0, (72) 

then T 2 (a;, 9) preserves V n „_i and the infinite flag of the polynomial spaces V m ,m, 
m = 0, 1, 2, . . . (cf. Cases J J. 5, 7.2.5, I5.^ } III. 4, 11.24 and II4.3). 
Case III. 2.1. If 

(2m — ?7)c +0 + nc +J + 2c + = , 
c+o - c+,7 = , 

mc , 5 — cy = 0, (73) 

at a certain integer m > 0, then T/2(x,0) preserves V n ,n-i and V m .rn-\- 
Case III. 2. 2. If 

c+o = c+,7 = c+ = , 

c+2 = cr=0 (74) 

then T 2 (x, 0) preserves 7-n,n-i ancl the infinite flag of polynomial spaces V m ,m-i> 
m = 0, 1, 2, . . . (cf. Cases 1.1.3, 1.2.3, 1.34, H.14, U.24, H4-3 and III. 1.2). This 
case corresponds to exactly-solvable operators Efc. 



In ]13j it has been shown that under a certain condition some quasi-exactly- 
solvable operators T2(x) in one real variable can preserve two polynomial spaces 
of different dimensions n and m (see Lemma 1.3). It has been shown that those 
quasi-exactly-solvable operators T 2 (x) can be represented through the generators 
of s? 2 (R) in a projectivized representation characterized either by the mark n or 
by the mark m. The above analysis shows that the quasi-exactly-solvable opera- 
tors T 2 (x, 9) in two variables (one real and one Grassmann) possess an extremely 
rich variety of internal properties. They are characterized by different structures 
of invariant sub-spaces. However, generically the quasi-exactly-solvable operators 
T 2 (x,#) can preserve either one, or two, or infinitely many polynomial spaces. For 
the latter, those operators become 'exactly-solvable' (see Cases 1.1.3, 1.2.3, 1.34, 
H.14, 11.24, II4.3 and III. 1.2) giving rise to eigenvalue problems (36) possess- 
ing infinite sequences of polynomial eigenfunctions. In general, for the two latter 
cases the interpretation of T 2 (x, 9) in terms of osp(2, 2) generators characterized by 
different marks does not exist, unlike the case of quasi-exactly-solvable operators in 
one real variable. The only exceptions are given by Case III. 2.1 and Case III. 2. 2. 



18 



ALEXANDER TURBINER. 

FIGURES. 3.0-3.III 



Newton diagrams describing invariant subspaces Vn,m of the second-order poly- 
nomials in the generators of osp(2, 2). The lower line corresponds to the part of 
the space of zero degree in 9 and the upper line of first degree in 9. The letters 
without brackets indicate the maximal degree of the polynomial in x. The letters 
in brackets indicate the maximal (or minimal) possible degree, if the degree can 
be varied. The thin line displays schematically the length of polynomial in x (the 
number of monomials). The thick line shows that the length of polynomial can 
not be more (or less) than that size. The dashed line means that the length of 
polynomial can take any size on this line. If the dashed line is unbounded, it means 
that the degree of the polynomial can be arbitrary up to infinity. The numbering 
of the figures I-III corresponds to the cases, which satisfy the conditions (45), (46) 
{Case I); (45), (47) (Case II) and (45), (48) {Case III). 
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Fig. 3.0. Basic subspace 
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Fig. 3.1.1 . Subspaces for the Case 1.1 
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Fig. 3.1.2 . Subspaces for the Case 1.2 
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Fig. 3.1.3 . Subspaces for the Case 1.3 
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Fig. 3.II.1 . Subspaces for the Case II. 1 
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Fig. 3.II.3 . Subspace for the Case IIS 
(!) m + l ( 2 ) m + 2 + k ( 2 ) m + 2 + k 



m m m 

(a) (b) (c) 

Fig. 3. II. 4 . Subspaces for the Case II. 4 
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Fig. 3.III.1 . Subspaces for the Case III.l 
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Fig. 3.III.2 . Subspaces for the Case III. 2 
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3.2. 2x2 matrix differential equations in x. It is well known that anti- 
commuting variables can be represented by matrices. In our case the matrix rep- 
resentation is as follows: substitute 9 and dg in the generators (31), (32) by the 
Pauli matrices a + and a~ , respectively, acting on two-component spinors. In fact, 
all main notations are preserved such as quasi-exactly-solvable and exactly-solvable 
operator, grading etc. 

In the explicit form the fermionic generators (32) in matrix representation are 
written as follows: 



Q = 



' a+- 




X(T+ 


, Q = 



xa d x — no~ 
-a~d x 



(75) 



The representation (75) assumes that in the spectral problem (36) an eigenfunction 
ip(x) is treated as a two-component spinor 



<p(x) 



'ipi(x)' 



In the matrix formalism, the polynomial space (29) has the form: 

Vn,m = 



JU j j • • • j Jli 

iXj j xL> •) • • • j *AJ 



(76) 



(77) 



where the terms of zero degree in 9 come in as the lower component and the terms 
of first degree in 9 come in as the upper component. The operator Tk(x, 9) becomes 
a 2 x 2 matrix differential operator T/,(x) having derivatives in x up to fc-th order. 
In order to distinguish the matrix operator in x from the operator in x, 9, we will 
denote the former as Tk(x). Finally, as a consequence of Theorem 3.1, we arrive at 
the eigenvalue problem for a 2 x 2 matrix quasi-exactly-solvable differential operator 
Tk(x), possessing in general (2n+l) polynomial solutions of the form V n n —\. This 
eigenvalue problem can be written in the form (cf. Eq.(37)) 



i—k 

H afe > 

i=0 



(x)ip x l) {x) = £lf(x) 



(78) 



where the notation (p x l ' means the ith order derivative with respect to x of each 
component of the spinor ip(x) (see Eq. (76)). The coefficient functions a.k t i(x) are 
given by 2 x 2 matrices and generically for the A:th order quasi-exactly-solvable 
operator their matrix elements are polynomials. Suppose that k > 0. Then the 
matrix elements are given by the following expressions 



*k,i( x ) 



k i 
n [k+i+l] 



B 



[k+i-l] 



k,i 
n [k+i 



(79) 



at k > 0, where the superscript in square brackets displays the order of the corre- 
sponding polynomial. 

It is easy to calculate the number of free parameters of a quasi-exactly-solvable 
operator Tk(x) 



par(T k (x))=4(k + l) 2 



(80) 



(cf. Eq.(40)). 
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For the case of exactly-solvable problems, the matrix elements (79) of the coef- 
ficient functions are modified 

«*,<(*) = f c ^i] J (81) 

where k > 0. An infinite family of orthogonal polynomials as eigenfunctions of 
Eq. (78), if they exist, will occur, if and only if the coefficient functions have the 
form (81). The number of free parameters of an exactly-solvable operator Efe(:r) 
and, correspondingly, the maximal number of free parameters of the 2x2 matrix 
orthogonal polynomials in one real variable, is equal to 

par(E k (x)) = 2fc(fc + 3) + 3 (82) 

(cf. Eq.(42)). 

The increase in the number of free parameters for the 2x2 matrix operators 
with respect to the case of the operators in x, 9 is connected to the occurrence of 
extra monomials of degree (fc + 1) in generators of osp(2,2) (see Eqs.(31), (32), 
(75)), leading to the kth order differential operators in x. 

Thus, the above formulas describe the coefficient functions of matrix differen- 
tial equations (78), which can possess polynomials in x as solutions, resolving the 
analogue of the generalized Bochner problem Eq. (0) for the case of 2 x 2 matrix 
differential equations in one real variable. 

Now let us take the quasi-exactly-solvable matrix operator T2(a;) and try to 
reduce Eq. (36) to the Schroedinger equation 

Id 2 
2dy 2 

where V(y) is a two-by-two hermitian matrix, by making a change of variable x i— > y 
and "gauge" transformation 



= E *(y) (83) 



* = XJip (84) 

where U is an arbitrary 2x2 matrix depending on the variable y. In order to get 
some "reasonable" Schroedinger equation one should fulfill two requirements: (i) 
the potential V(y) must be hermitian and (ii) the eigenfunctions ^(y) must belong 
to a certain Hubert space. 

Unlike the case of quasi-exactly-solvablc differential operators in one real variable 
(see H ) , this problem has no complete solution so far. Therefore it seems instructive 
to display a particular example fll4| . 

Consider the quasi-exactly-solvable operator 

T 2 = -2T°T- + 2T-J - i/?T°Qi + 

aT° - (2n + 1)T~ - y (3n + 1)Q J + % -a{5Q 2 - i0Q x , (85) 

where a and (3 are parameters. Upon introducing a new variable y = x 2 and after 
straightforward calculations one finds the following expression for the matrix U in 
Eq. (84) 

U = cxp( — + —£-<ri) (86) 



22 



ALEXANDER TURBINER 



and for the potential V in Eq. (83) 

V(y) = \{a 2 - (3 2 )y 2 + a 2 [-(n + \)(J + f y 2 - \ tan if] cos ^ + 

^3h("+ i )/3+^2; 2 - I cot^|-]sin^- (87) 

It is easy to see that the potential V is hermitian; (2n + 1) eigenfunctions have the 
form of polynomials multiplied by the exponential factor U and they are obviously 
normalizable. 

4. Polynomials in two real variables 

For further consideration it is convenient to illustrate the space of polynomials 
of finite degree in several variables through Newton diagrams. In order to do this, 
let us introduce a d-dimensional integer lattice in R d and put into correspondence 
to each node with the coordinates (fci, fe, • ■ ■ kd) the monomial xi kl X2 k2 ■ ■ ■ Xd kd ■ 

4.1. Polynomials of the first type. Now let us describe the projectivized repre- 
sentation of the algebra sl^ (R) in the differential operators of the first order acting 
on functions of two real variables. It is easy to show that the generators have the 
form ft| 

J3 = y 2 d y + xyd x - ny , j\ = x 2 d x + xyd y - nx , 
J 2 = -yd x , J 2 = -d x , J'f = -d y , Jf = -xd v , 

Jd = yd y + 2xd x - n , Jd = 2yd y + xd x - n , (88) 

where x, y are the real variables and n is a real number. If n is a non-negative 
integer, the representation becomes finite-dimensional of dimension (l+n)(l + n/2). 
The invariant sub-space has a polynomial basis and is presented as a space of all 
polynomials of the following type 

7>« = (1; x, y- x 2 , xy, y 2 ; . . . ; x n , x^y, . . . , xy"^ 1 , y n ) (89) 

or, graphically, the space (89) is given by the Newton diagram of Figure 4.1. 




Fig. 4.1 . Graphical representation (Newton diagram) of the space Vn 

(see (89)). 

Definition 4.1. Let us name a linear differential operator of the fc-th order a 
quasi-exactly-solvable of the first type,T k {x, y), if it preserves the space Vn^ ■ 
Correspondingly, the operator (x, y) € [x, y), which preserves the infinite 
flag Vq 1 ^ C C C ■ • • C Vn^ C . . . of spaces of all polynomials of the type 
(89) , is named an exactly-solvable of the first type. 
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Lemma 4.1. Consider the space "Pn ^ ■ 

(i) Suppose n > (k — 1). Any quasi-exactly-solvable operator of the first type 
T^\x,y), can be represented by a k-th degree polynomial of the generators (88). 
If n < (k — 1), the part of the quasi-exactly-solvable operator T^\x,y) of the first 
type containing derivatives in x, y up to the order n can be represented by a n-th 
degree polynomial in the generators (88). 

(ii) Conversely, any polynomial in (88) is a quasi-exactly solvable operator of 
the first type. 

(Hi) Among quasi-exactly-solvable operators of the first type there exist exactly- 
solvable operators of the first type E^^y) C (x, y). 

Comment 5. The meaning of the lemma is the following: (x, y) at k < n+l is 
simply an element of the universal enveloping algebra U s [ 3 ^ of the algebra s^R) 
in realization (88) . If k > n+l, then (x, y) is represented as a polynomial of nth 
degree in (88) plus B dx , n - m +i dym , where m = 0, 1, . . . (n + 1) and B is any linear 
differential operator of order not higher than (k — n — 1). 

Let us introduce the grading of the generators (88) in the following way. The 
generators are characterized by the two-dimensional grading vectors 

deg{Jl) = (0,+l) , deg(Jl) = (+1,0) , 

deg(Jl) = (-1.+1) , deg{4) = (+1,-1) , 

deg(J d ) = (0,0) , deg{J d ) = (0,0) , 

deg{Jl) = (0,-1) , deg(Jl) = (-1,0). (90) 

It is apparent that, the grading vector of a monomial in the generators (88) can be 
defined by the grading vectors of the generators by the rule 

d~eg(T) = ^[(Ja 1 )" 1 ^^)" 1 ^^)" 2 ^^)" 3 ^^)"^^)"^^)" 3 ^^)" 21 ] = 
("13 + "32 - "23 - "31 , "12 + "23 - "32 - "21 ) = (deg x (T) , deg y (T)) 

(91) 

Here the n's can be arbitrary non-negative integers. 

Definition 4.2. Let us name the grading of a monomial T in generators (75) the 
number 

deg{T) = deg x {T) + deg y {T) . 

We will say that a monomial T possesses positive grading if this number is positive. 
If this number is zero, then a monomial has zero grading. The notion of grading 
allows one to classify the operators T^\x,y) in a Lie-algebraic sense. 

Lemma 4.2. A quasi-exactly-solvable operator C C^ s i 3 (R) has no terms of 

positive grading if and only if it is an exactly-solvable operator of the first type. 

It is worth noting that among exactly-solvable operators there exists a certain 
important class of degenerate operators, which preserve an infinite flag of spaces of 
all homogeneous polynomials 

V£> = (x m ,x m - 1 y,...,xy m - 1 ,y m ) (92) 

(represented by a horizontal line in Fig. 4.1). 
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Lemma 4.3. A linear differential operator Tfe(x, y) preserves the infinite flag C 
C C • • • C VrP C . . . 0/ spaces of all polynomials of the type (89), if 
and only if it is an exactly- solvable operator having terms of zero grading only. 
Any operators of such a type can be represented as a polynomial in the generators 
Ji,Ji,Jd,Jd ( see (88)), which form the algebra so 3 © R . If such an operator 
contains only terms with zero grading vectors, this operator preserves any space of 
polynomials. 

Theorem 4.1. Let n be a non-negative integer. In general, the eigenvalue problem 
for a linear symmetric differential operator in two real variables Tk(x,y): 

T k {x,y)tp(x,y) = eip(x,y) (93) 

has (n + l)(n/2 + 1) eigenf unctions in the form of a polynomial in variables x,y 
belonging to the space (89), if and only if Tk(x,y) is a symmetric quasi- exactly- 
solvable operator of the first type. The problem (93) has an infinite sequence of 
eigenf unctions in the form of polynomials of the form (89), if and only if the oper- 
ator is a symmetric exactly- solvable operator. 

This theorem gives a general classification of differential equations 

m=0 i=0 ' y 

having at least one eigenfunction in the form of a polynomial in x, y of the type (89). 
In general, the coefficient functions a^ l _ i (x,y) have quite cumbersome functional 
structure and we do not display them here (below we will give their explicit form 
for T^\x, y)). They are polynomials in x, y of order (k + m) and always contain 
a general inhomogeneous polynomial of order to as a part. The explicit expres- 
sions for those polynomials are obtained by substituting (88) into a general, fcth 
order polynomial element of the universal enveloping algebra C/ s i 3 (R) of the algebra 

sZ 3 (R.) . Thus, the coefficients in the polynomials a\"^_ i (x,y) can be expressed 
through the coefficients of the k-th order polynomial clement of the universal en- 
veloping algebra C/ s ; 3 (r). The number of free parameters of the polynomial solu- 
tions is defined by the number of parameters characterizing a general, fcth order 
polynomial element of the universal enveloping algebra C/ s i 3 (R). In counting free 
parameters some relations between generators should be taken into account, specif- 
ically for the given representation (88) 

J^^d ~ 2J 2 Jd — SJgjf = n ^2 
J 3 Jd 2J^J d 3J^J^ — ^^3 

Jfj d + JfJd - ZJ^Jl = (n + 3)j| 
JlJ d + J\Jd - 3J\j\ = {n + 3)Jl 
3(J 2 1 J 1 2 + J\JI + J|J 3 2 ) + J d J d +J d J d - JdJd = 3J d + 3n + n 2 
2J d J d + 2J d J d - 5J d J d + 9J 2 3 J 3 2 = (n + 6)J d + (n - 3) J d + n 2 + 3n 
JdJd - JdJd + 3{JlJl - J}Jl) = (n + 3)(J d - J d ) 
J d j\ — 2J d j\ — 3J3J1 = nj\ 



J d j\ — 2J d j\ — 3J2J 2 = njf 



(95) 
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between quadratic expressions in generators (and the ideals generated by them).^ 
For the case of exactly-solvable problems, the coefficient functions a!f^_ i {x 1 y) take 
the form 

p+t/<m 

a tm-i( x >y) = a ^m- l , P , q x p y q (96) 

p,q=0 

with arbitrary coefficients. 

Now let us proceed to the case of the second-order differential equations. The 
second-order polynomial in the generators (88) can be represented as 

T% = C a /3 l7 sJpJ] + C a pJp + c , 

where we imply summation over all repeating indices; a,/3, 7,(5 correspond to the 
indices of operators in (88) and for the Cartan generators we suppose both indices 
simulate d or d, all c's are set to be real numbers. After substitution the expressions 
into T 2 the explicit form of the quasi-exactly-solvable operator is given by 

[x 2 P^(x,y) + xP%(x,y) + P™(x,y)}-^ + 

d 2 



[xyP%Ux,y) + P*\{x,y) + P* y (x,y)] 



dxdy 

[y 2 P yy 2 (x,y) + yP^{x,y) + pwfay)]^ + 
[xPZ <2 (x,y) + P^(x,y) + p- ( x ,y)]^- + 
[yPl 2 (x,y) + P^y) + 

[P° 2 (x,y) + P? tl (x,y) + P§ t0 (x,y)] (97) 

where -P^ m (^,y) anc ^ Pkm( x >y) are homogeneous and inhomogeneous polynomials 
of order k, respectively, the index m numerates them, the superscript c characterizes 
the order of derivative, which this coefficient function corresponds to. For the case 
of the second-order-exactly-solvable operator 'E 2 I \x, y), the structure of the coeffi- 
cient function is similar to (97), except for the fact that all tilde-less polynomials 
disappear. The number of free parameters is equal to 

par(TP(x,y))=36. 

while for the case of an exactly-solvable operator (an infinite sequence of polynomial 
eigenfunctions in the problem (93)) 

par(E { 2 I \x,y)) = 25 . 

An important particular case is when the quasi-exactly-solvable operator (x, y) 
possesses two invariant sub-spaces of the type (89). This situation is described by 
the following lemma: 



7 It was shown that there exist 56 cubic relations between the generators (88). However, it 
turned out that all of them (including standard cubic Casimir operator in given representation 
(88)) are functionally-dependent on quadratic relations (95) (Personal communication by G. Post). 
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Lemma 4.4. Suppose T 2 ( x iV) has no terms of grading 2 

Cl2,12 = Ci 3l i3 = Ci2,13 = . (98) 

If there exist some coefficient c 's and a non-negative integer N such that the con- 
ditions 

C12 = (n - N - m)ci2,d + (n - 2N + m)c 12 j + (N - m)c 12 ^2 , 

C13 = (n - N - m)ci 3t d + {n - 2N + m)c 13 j + toci 3;2 3 , (99) 

are fulfilled at all m = 0,1,2, ... , N , then the operator { x i v) preserves both V n 
andVyi, and par(T 2 (x,y)) — 31. 

Now let us proceed to the important item: under what conditions on the coeffi- 
cients in (97) the second-order-quasi-exactly-solvable operators can be reduced to 
a form of the Schroedinger operator after some gauge transformation (conjugation) 

f{x,y)e t ^T 2 {x,y)^e- t ^ = -A g + V(x,y,n) (100) 

where f,t,V are some functions in R 2 , and A g is the Laplace-Beltrami operator 
with some metric tensor g^ u \ (j,i/= 1,2 0, This is a difficult problem which has 
no complete solution yet. In |l4|, |l6) a few multi-parametric examples have been 
constructed. 

In general, if one can perform the transformation (100), we arrive at potentials 
V(x, y, n) containing explicitly a dependence on the parameter n and hence the 
dimension of the initial invariant subspace (89). Finally, we end up at quasi- exactly- 
solvable Schroedinger equations with a certain number of eigenfunctions known 
algebraically. Taking an exactly-solvable operator Eg (x, y) instead of the quasi- 
exactly-solvable one and performing the transformation (100), a certain exactly- 
solvable Schroedinger operator emerges. Evidently, they will have no dependence 
on the parameter n and correspondently an infinite sequence of eigenfunctions of 
the Schroedinger equation has a form of polynomials. 

Also it is worth noting that there exists an important particular class of quasi- 
exactly-solvable operators, lying in a certain sense in between quasi-exactly-solvable 
and exactly-solvable operators. The algebra s?3(R) in realization (88) contains the 
algebra so 3 (R) as a sub-algebra in such a way that 

J 1 = (1 + y 2 )d y + xyd x — ny , J 2 = (1 + x 2 )d x + xyd y — nx , 

J 3 = xd y - yd x (101) 

which is not graded anymore. If the parameter n is a non-negative integer (and co- 
incides with that in (88)), the same finite-dimensional invariant sub-space Vn^ (see 
(89) and Fig. 4.1) as in the original sis(R) occurs. For this case a 
corresponding finite-dimensional representation is reducible and unitary. It has 
been proven JlJ, 0, that any symmetric bilinear combination of generators (101), 
T 2 = c a /3J a J l3 , where c a p = cp a are real numbers, can be reduced to a form of the 



Our further consideration will be restricted the case / = 1 only. 
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Laplace-Beltrami operator plus a scalar function making a certain gauge transfor- 
mation (conjugation). [] In general, the metric tensor g^ v is not degenerate ^| and 
the obtained potential V(x,y) is given by a rational function of the coordinates. 
Moreover, the potential V(x, y) has no dependence on the label nl All n-dependent 
terms contain no x, y- dependence and can be included into redefinition of the ref- 
erence point for eigenvalues. At first sight, those quasi-exactly-solvable operators 
look as if they are exactly-solvable, but is not so. As a consequence of the fact 
that the quadratic Casimir operator for sc>3(R) in the realization (101) is non- 

(e) . (e) 

trivial and commutes with T 2 i the functional space of T 2 is subdivided into the 
finite-dimensional blocks of increasing sizes, corresponding to the irreducible rep- 
resentations of S03(R) (unlike truly exactly-solvable operators, for which the sizes 
of blocks are fixed and equal to one). We name such operators exactly- solvable of 
the second type. 

Another interesting property appears, if the matrix c a g has one vanishing eigen- 
value: a certain mysterious relation holds H] (see also Jl7], Q ) 

V(x,y,{c}) = ^R(x,y,{c}), (102) 

where R(x,y,{c}) is the scalar curvature calculated through the metric tensor 
attached in the Laplace-Beltrami operator. This can imply that the correspond- 
ing Schroedinger operator has the form of a purely geometrical object! The real 
meaning of this fact is not understood so far. 

4.2. Polynomials of the second type. In jll| we studied the quasi-exactly- 
solvable operators in one real variable. It turned out that the solution to this 
problem was found using a connection with the projectivized representation of the 
algebra s^(R)- As a natural step in developing the original idea, let us consider 
the projectivized representation of the direct sum of two algebras 5/2 (R)- 

The algebra s^(R) © s^(R) taken in projectivized representation acts on func- 
tions of two real variables. The generators have the form (see (2)) 

J x = x 2 d x - nx , Jy = y 2 d y - my , 

tO o n j0 o m 

j x — xo x — — , J y — yOy — , 

J x =d x , Jy = d y , (103) 

where x, y are the real variables and n, m are numbers. If n, m are non-negative inte- 
gers, there exists the finite-dimensional representation of dimension 
(n + l)(m + 1). Evidently, the invariant sub-space has a polynomial basis and 
is presented as a space of all polynomials described by the Newton diagram of Fig. 
4.2. We denote this space as Vn,m- 



9 In fact, in |17| has been proven a more general theorem, that any symmetric bilinear com- 
bination of generators of unitary representation of the semisimple Lie algebra, realized in the 
first-order differential operators, can be reduced to a form of the Laplace-Beltrami operator with 
some metric tensor plus a scalar function making a certain gauge transformation (conjugation). 

10 Degeneracy of g M „ occurs, for example, if c a p has two vanishing eigenvalues. 




Fig. 4.2. Graphical representation (Newton diagram) of the space Vn.m ■ 

Definition 4.3. Let us name a linear differential operator of the (fc,p)th order, 
containing derivatives in x and y up to fcth and pth orders, respectively^^, a 
quasi-exactly-solvable of the second type, (x,y), if it preserves the Space T^n.rn- 

Correspondingly, the operator Ejj.^ (x, y) S (x,y), which preserves either the 
infinite flag V^JJ C V { "2 C C ■ ■ • C vi'l C ... , or the infinite flag 

Pnfi c ^nfi c 7*71,2 C • ■ ■ C PrwL C . . . of spaces of all polynomials, is named 
an exactly-solvable of the type 2 X or 2 V , respectively. 

Lemma 4.5. Consider the space Vn.m- 

(i) Suppose n > (k—1) and m > (p — 1). Any quasi-exactly-solvable operator of 
the second type ^ (x,y), can be represented by a (k,p)th degree polynomial of the 
generators (103). If n < (k—1) and/or m < (p — 1), the part of the quasi-exactly- 
solvable operator Tffp (x,y) of the second type containing derivatives in x,y up to 
order n,m, respectively, can be represented by an (n,m)th degree polynomial in the 
generators (103). 

(ii) Conversely, any polynomial in (103) is a quasi-exactly solvable operator of 
the second type. 

(Hi) Among quasi- exactly- solvable operators of the second type there exist exactly- 
solvable operators of the second type (x, y) C (x,y). 

Similarly, as it has been done for the algebra s^(R), one can introduce the notion 
of grading: 

deg(J+) = (+1,0) , deg(J+) = (0,+l) , 
deg(J° x ) = (0,0) , deg(J Q y ) = (0,0) , 

deg(J-) = (-1, 0) , deg(J~) = (0, -1). (104) 

It is apparent that, the grading vector of a monomial in the generators (90) can be 
defined by the grading vectors (104) of the generators by the rule 

deg(T) = deg[(J+)^(J°)^(J^) n '-(J+)^(J°)^(J-) n '-] = 

(n x+ - n x - , n y+ - n y -) = (deg x (T) , deg y (T)) (105) 
Here the n's can be arbitrary non-negative integers. 



So a leading derivative has a form k p . Also we will use a notation through this section 
Tpf{ x , y) implying that in general all derivatives of the order N are presented. 
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Definition 4.4. Let us name the .T-grading, y-grading and grading of a mono- 
mial T in generators (103) the numbers deg x (T), deg y (T) and deg(T) — deg x (T) + 
deg y (T), respectively. We say that a monomial T possesses positive x- grading 
(y-grading, grading), if the number deg x (T) (deg y (T), deg(T)) is positive. If 
deg x (T)(deg y (T),deg(T)) = 0, then a monomial has zero x-grading (y-grading, 
grading). 

The notion of grading allows one to classify the operators T^*\x,y) in the 
Lie-algebraic sense. 

Lemma 4.6. The quasi-exactly-solvable operator ^ (x , y) preserves the infinite 

fl a 9 T'ofm c 'P'l.m c T^fm C ■ ■ ■ C Vnjm C . . . of spaces of all the polynomials, if 
and only if it is an exactly-solvable operator of type 2 X having no terms of positive 
x-grading, deg x (T) > 0. 

The quasi-exactly-solvable operator T^\x,y) preserves the infinite flag V^fi c 

VnP C VnP C • • • C Vn,m C . . . of all spaces of the polynomials, if and only if 
it is an exactly- solvable operator of type 2 y having no terms of positive y-grading, 
deg y {T) > 0. 

If a quasi-exactly-solvable operator of the second type contains no terms of 
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positive grading, this operator preserves the infinite flag V^ I] C v[ I] C C • • • C 
Vn^ C . . . of spaces of all the polynomials of the type (76) and is attached to the 
exactly-solvable operator of the first type. 

Theorem 4.2. Let n,m be non-negative integers. In general, the eigenvalue prob- 
lem (93) for a linear symmetric differential operator in two real variables Tk,p(x, y) 
has (n + l)(m + 1) eigenf unctions in the form of a polynomial in variables x, y be- 
longing to the space Vn,m , if and only if T^ p {x, y) is a quasi- exactly- solvable 
symmetric operator of the second type. The problem (93) has an infinite sequence 
of eigenf unctions in the form of polynomials belonging the space Vn,m at fixed to 
(n), if and only if the operator is an exactly-solvable symmetric operator of type 2 X 
fty). 

This theorem gives a general classification of differential equations (94), hav- 
ing at least one eigenfunction in the form of a polynomial in x,y of type Vn.rl- 
In general, the coefficient functions a["^_ i (x,y) in (94) have a quite cumbersome 
functional structure and we do not display them here (below we will give their 
explicit form for Tj ( x iV))- They are polynomials in x,y of the order (fc + to). 
The explicit expressions for those polynomials are obtained by substituting (103) 
into a general, fcth order polynomial element of the universal enveloping algebra 
U s i 2 (r)®si 2 (r) of the algebra s?2(R) © s?2(R) ■ Thus, the coefficients in the polyno- 
mials af^_ i (x,y) can be expressed through the coefficients of the k-th order poly- 
nomial element of the universal enveloping algebra C/ s ; 2 (r)© s / 2 (r) ■ The number of 
free parameters of the polynomial solutions is defined by the number of parameters 
characterizing a general, kth order polynomial clement of the universal enveloping 
algebra C/ s ; 2 (r)© s / 2 (r). In counting free parameters some relations between gener- 
ators should be taken into account, specifically for the given representation (103) 

I / \ 

Jx Jx ~ Jx"x "f" x = — ~2\~2 ' 

J+J- - J°J° + J° = -^(^ + 1) (106) 

between quadratic expressions in generators (and the ideals generated by them) . ^ 
Now let us proceed to the case of the second-order differential equations. The 
second-order polynomial in the generators (103) can be represented as 



To — r xx T a T 13 -4- r xy T a T@ 4- r vv T a T 13 -L r x T a J- r y T a 4- r 
±2 — ^aj3 u x J x < u a0" J x u y ' aj3' J y ° y ' ^a°x ' ^a J y ' L 



(107) 



where we imply summation over all repeating indices and a, (3 — ±,0; all c's are 
set to be real numbers. Taking (106) in account, it is easy to show that T2 is 
characterized by 26 free parameters. Substituting (103) into (107), we obtain the 
explicit form of the second-order-quasi-exactly-solvable operator 

T { 2 H \x,y) = 



For this case they correspond to the quadratic Casimir operators. 
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Ox dy (1Q8) 

where P£ m (x,y) anc ^ Pkm( x ^y) are homogeneous and inhomogeneous polynomials 
of order k, respectively, the index m marks them, and the superscript c characterizes 
the order of the derivative corresponding to this coefficient function. For the case 
of the second-order-exactly-solvable operator of type 2 X 

Qm*)-^ + WZ(y) + + + 

Qf,iW + yQU x )}w- + QIM + *Q?,o(«)]#- + Q%(y) 

Ox dy (10g) 

For the case of the second-order-exactly-solvable operator of type 2 y the functional 
form is similar to (109) with the interchange x *-* y. 

It is easy to show that the numbers of free parameters are 

par(Tf 7) (a;,y)) = 26 . 

par(E { 2 H) (x,y)) = 20 . 

As in the previous cases, there is an important particular case of quasi-exactly- 
solvable operators of the second order, where they possess two invariant sub-spaces. 

Lemma 4.7. Suppose in (107) there are no terms of ' x- grading 2 

= 0. (110) 



If there exist some coefficient c 's and a non-negative integer N such that the 
conditions 

c x A = o , 



, 



? + = (n/2 - N)c%l + (m/2 - k)c x + v , (111) 

are fulfilled at all k = 0,1,2, ... , to, then the operator T^ T \x, y) preserves both 
V% m andPfl^, and par{T ( £ I \x,y)) = 22. 

Generically, the question of the reduction of quasi-exactly-solvable operators of 
the second type T2 (a;, y) to the form of the Schroedinger operator is still open. 
Recently, several multi-parametrical families of those Schroedinger operators were 
constructed in |l4|, |l6). As well as the case of the first type of quasi-exactly- 
solvable operators, corresponding Schroedinger operators contain in general the 
Laplace-Beltrami operator characterizing a non-flat space metric tensor. 

The above analysis of linear differential operators preserving the space Vn,m 
can be naturally extended to the case of linear finite-difference operators defined 
through the Jackson symbol D (see Section 2) with the action: 

Df(x) = fi f- f ^ + f(ax)D 
(1 - q)x 
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instead of a continuous derivative, where f(x) is a real function and q is a number. 
All above-described results hold with replacement of the algebra s/2 (R) © sl 2 (R) 
by the quantum algebra sl2(R) q © s?2(R)<j in 'projectivized' representation (21) 
(see [|| |3)). Obviously, in the limit q — > 1 all results coincide to the results of the 
present Section. 



4.3. Polynomials of the third type. The third case, which we shall discuss here, 
corresponds to the family of Lie algebras gl% (R) tx R r+1 (semidirect sum of gl^jR) 
with a (r + l)-dimensional abelian ideal; Case 24 in the classification given in |19[| ). 
This family of Lie algebras, depending on an integer r > 0, can be realized in terms 
of the first-order differential operators 

J 1 = d x , 
J 2 = xd x - I , J 3 = ydy - — , 
J 4 = a; 2 d x + rxydy — nx , 



J 



5+i 



= x l d y , i = 0, 1, 



(112) 



where a;,y are real variables and n is a real number. |^ If n is a non-negative 
integer, the representation becomes finite-dimensional. The invariant sub-space 
has a polynomial basis and is presented as a space of all polynomials of the form 



i-\-rj<n 
i,J>0 



(113) 



or, graphically, (113) is given by the Newton diagram of Figure 4.3. The general 
formula for the dimension of the corresponding finite-dimensional representation 
(113) is given by 



where for small r 



diraV, 



(in) _ [n 2 + (r + 2)n + a ri n] 
2r 



Oil- 



(114) 



0-2, 



Ol3,; 



3 at odd n 

4 at even n 



4 at (n+l)=0 (mod 3) 
6 at other n 



5 at (n+l)=0 (mod 4) 
ot-A.n — \ 8 at other n 

9 at (n+3)=0 (mod 4) 



13 with minor modifications 

14 It is worth noting that at r = 1 the algebra {gl-2(R) IX R 2 } C s/3(R). Thus, this case is 
reduced to one about the first type polynomials (see Section 4.1). Hereafter, we include r = 1 
into consideration just for the sake of completeness. 
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Fig. 4.3 . Graphical representation (Newton diagram) of the space 

,n 

(see (113)). 

Definition 4.5. Let us name a linear differential operator of the fcth order a quasi- 
exactly-solvable of the r -third type, T < j^' III \x,y), if it preserves the space V^n 1 ^ ■ 
Correspondingly, the operator 'E^' III \x, y) € T^' III \x,y), which preserves the 
infinite flag P<Jf 7) C P$." n C P^" z) C • • • C P$ z) C ... of spaces of all polyno- 
mials of the type (113) , is named an exactly- solvable of the r-third type. 

Lemma 4.8. Take the space V^n^ ■ 

(i) Suppose n > (fc — 1). Any quasi-exactly-solvable operator r T < ^' III \x,y) can 
be represented by a kth degree polynomial of the generators (112). If n < (k — 
1), the part of the quasi- exactly-solvable operator T < j^' III \x,y) of the r-third type 
containing derivatives in x, y up to order n can be represented by a n-th degree 
polynomial in the generators (112). 

(ii) Conversely, any polynomial in (112) is a quasi-exactly solvable operator of 
the r-third type. 

(Hi) Among quasi-exactly-solvable operators of the r-third type there exist exactly- 
solvable operators of the r-third type E,'^' III \x,y) C T < ^' III \x,y). 

One can introduce the grading of the generators (112) in an analogous way as has 
been done before for the cases of the algebra s^(R) (see (90)) and s^R) © s^(R) 
(see (104)). All generators are characterized by the two-dimensional grading vectors 

cie ff (J 1 ) = (-l,0) , 
deg(J 2 ) = (0,0) , cfe 5 (J 3 ) = (0,0) , 
deg(J 4 ) = (1,0) , 

deg(J 5 ) = (0, -1) , deg(J 6 ) = (1, -1) , . . . , deg( J 5+r ) = (r, -1) 

(115) 

Similarly as before, the grading vector of a monomial in the generators can be 
defined through the grading vectors of the generators (112) (cf. (90), (105)). 

Definition 4.6. Let us name the grading of a monomial T in generators (112) the 
number 

deg(T) = deg x (T) + r deg v {T) 
(cf. the case of s^(R)). We say that a monomial T possesses positive grading if 
deg(T) is positive. If it is zero, then a monomial has zero grading. 

The notion of grading allows us to classify the operators Tk(x,y) in a 
Lie-algebraic sense. 
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Lemma 4.9. A quasi-exactly-solvable operator T^' 111 ^ c f/ g j 2 (RWRr+i has no 
terms of positive grading iff it is an exactly-solvable operator. 

Theorem 4.3. Let n and (r — 1) be non-negative integers. In general, the eigen- 
value problem (93) for a linear symmetric differential operator in two real variables 
Tk(x,y) has a certain number of eigenf unctions in the form of a polynomial in 
variables x,y belonging to the space (113), if and only ifT/,(x,y) is a quasi- exactly- 
solvable, symmetric operator of the r -third type. The problem (93) has an infinite 
sequence of eigenf unctions in the form of polynomials belonging to (113), if and 
only if the operator is an exactly-solvable, symmetric operator of the r -third type. 

This theorem gives a general classification of differential equations (94), having 
at least one eigenfunction in the form of a polynomial in x,y of the type P^n^ ■ 
In general, the coefficient functions a\"^_ i (x,y) in (94) are polynomials in x,y and 
have a quite cumbersome functional structure and we do not display them here 
(below we will give their explicit form for Tf, r,/// ' (a;, y)). The explicit expressions 
for those polynomials are obtained by substituting (112) into a general, fcth order 
polynomial clement of the universal enveloping algebra ^ s i 2 (r)kR'-+i of the algebra 

gh(R) x R r+1 . Thus, the coefficients in the polynomials a\"^_ i (x,y) can be ex- 
pressed through the coefficients of the fcth order polynomial element of the universal 
enveloping algebra £/ ff i 2 (R)xR r- + 1 - The number of free parameters of the polynomial 
solutions is defined by the number of parameters characterizing a general, kth order 
polynomial element of the universal enveloping algebra. In counting free param- 
eters some relations between generators should be taken into account, specifically 
for a given representation (112) 

J 2 J 5 - J 1 J 6 + £ J 5 = , 



3 

J 1 J 4 - J 2 J 2 - rJ 2 J 3 — J 2 — r(^ + l)J 3 = -|(| + 1) 



J 2 J 6+t + rJ 3 J 6 +* - J 4 J 5 +* - (- + 1) J 6+t = , 

3 



(116) 



at i = 0,1,2,..., (r-1) , (117) 

J 1 J 7+i - J 2 J 6 +* - (- + 1) J 6+l = , 
3 

at i = 0,l,2,...,(r-2) , (118) 

j5j7+i _ . . . _ j5+k j7+i—k 

at k = 0,1,2,... and 2k < (2 + i) , and i = 0, 1, 2, . . . , (2r - 2) 

(119) 

between quadratic expressions in generators (and the ideals generated by them). 

Now let us proceed to the case of the second-order differential equations. Taking 
the relations (116)-(119) in account, one can find the number of free parameters 

par{T^ nl \x,y)) = 5(r + 4) (120) 
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and obtain the explicit form of the second-order-quasi-exactly-solvable operator 

[PZA x ) + ftuW + y^.oW]^ + p U x ) ( 121 ) 

where Pu m {x, y) are inhomogeneous polynomials of order k, the index m numerates 
them, and the superscript c characterizes the order of the derivative correspond- 
ing to this coefficient function. For the case of the second-order-exactly-solvable 
operator of r-third type 

ET ,m) (*,3/) = 

p »)^ + + y p ^^d- y + + ^!oW]]^ + 

[P£o(*) + Ri(^) + ^o,oW]^ + n°oW (122) 
In this case the number of free parameters is equal to 

par(V ( { JH \x,y)) = 5(r + 3) (123) 

As in the previous cases, there is an important particular case of r-third type quasi- 
exactly-solvable operators, where they possess two invariant sub-spaces. 

Lemma 4.10. Suppose E^'' III \x, y) has no terms of grading 2 

c 44 = . (124) 

If there exist some coefficient c 's and a non-negative integer N such that the con- 
ditions 

C4,5+r = , 

N N 
c 4 = (— - m)c 2 4 + (y - fc)c 34 , (125) 

are fulfilled at all m,k = 0,1,2,... such that m + rk = N, then the operator 
T l ^' III \x, y) preserves both V* 1 ^ and V*™, and par( r £^' III \x, y)) — 5r + 17. 

Generically, the question of the reduction of the quasi-exactly-solvable operator 
Tn"' III \x 1 y) to the form of the Schroedinger operator is still open. Initially in 



|16| a few multi-parametrical families of those Schroedinger operators have been 
constructed. Similar to the case of quasi-exactly-solvable operators of the first 
and the second types, corresponding Schroedinger operators contain in general the 
Laplace-Bcltrami operator with non-flat space metric tensor. 
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4.4. Discussion. Through all the above analysis, one very crucial requirement 
played a role: we considered the problem in general position, assuming the spaces of 
all polynomials contain polynomials with all possible real (complex) coefficients. If 
this requirement is not fulfilled, then the above connection of the finite-dimensional 
space of all polynomials and the finite-dimensional representation of the Lie algebra 
is lost and degenerate cases occur. This demands a separate investigation. 

So, we described three types of finite-dimensional spaces (see Fig. 4.1-3) with a 
basis in polynomials of two real variables, which can be preserved by linear differen- 
tial operators. The natural question emerges: is it possible to find linear differential 
operators having no representation through generators of a Lie algebra, which pre- 
serve the finite-dimensional space of polynomials, other than those shown in Fig. 
4.1,2,3? The answer is given by the following conjecture pC| . 

Conjecture 1. // a linear differential operator Tk{x,y) does preserve a finite- 
dimensional space of all polynomials (presented by a convex Newton diagram) other 
than Vn^ , Vn,m > "Pr*n , then this operator preserves either one of those spaces as 
well, or a certain infinite flag of finite- dimensional spaces of all polynomials. 

The general question about reducibility of a second-order polynomial in gener- 
ators of the algebras sZ 3 (R), sZ 2 (R) © sZ 2 (R) and {gh(B.) x R r+1 }, r = 2, 3, 4, . . . 
to the form of the Laplace-Beltrami operator plus a potential is still open. For all 
cases where this reducibility has been performed, a metric tensor appeared which 
corresponded to non-flat space. All attempts to find the quasi-exactly-solvable 
or exactly-solvable problems in R 2 corresponding to the flat space did lead to a 
problem with separability of variables. We believe that the following conjecture 
holds. 

Conjecture 2. In R 2 there exist no quasi-exactly-solvable or exactly-solvable 
problems containing the Laplace-Beltrami operator with flat-space metric tensor, 
which are characterized by non-separable variables. 

Also the next question, under what conditions the cigenfunctions obtained belong 
to some Z/2-space, is lacking a solution. The analysis analogous to that which has 
been performed in |l6[ ] for one-dimensional quasi-exactly-solvable problems should 
be done. 

5. General considerations 
As the main conclusion one can formulate as the following theorem. 

(Main) Theorem. Consider a Lie algebra g of the differential operators of 
the first order, which possesses a finite- dimensional irreducible representation V . 
Any linear differential operator acting in V (having V as the invariant subspace) 
can be represented by a polynomial in generators of the algebra g plus an operator 
annihilating V . 



As a probable extension of this theorem, we assume that the following conjecture 
holds 
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Conjecture 3. If a linear differential operator T acting on functions in R k 
does possess a single finite- dimensional invariant subspace with polynomial basis 
(presented by a convex Newton diagram), this finite- dimensional space coincides 
with a certain finite- dimensional representation of some Lie algebra, realized by 
differential operators of the first order. 

The interesting question is how to describe finite modules of smooth functions 
in R k which can serve as invariant sub-spaces of linear operators. Evidently, if 
we could find some realization of a Lie algebra g of differential operators, possess- 
ing an irreducible finite module of smooth functions, we can immediately consider 
the direct sum of several species of g acting on the functions given at the corre- 
sponding direct products of spaces, where the original algebra g acts in each of 
them. As an example, this procedure is presented in Section 4.2: the direct sum 
of two s^(R) generates the linear differential operators acting on the functions 
at R 2 having a rectangular Vnjn as the invariant sub-space. Taking the direct 
sum of k species of s^R) acting at R k , we arrive to the operators possessing a 
/c-dimensional parallelopiped as the invariant sub-space. In this case the algebra 
sl2 (R) plays the role of a "primary" algebra giving rise to a multidimensional con- 
vex geometrical figure (Newton diagram) as the invariant sub-space of some linear 
differential operators. In the space R 1 there is only one such "irreducible" convex 
Newton diagram - a finite interval, in other words, a space of all polynomials in x 
of degree not higher than a certain positive integer. Thus, in R 2 the rectangular 
Newton diagram is "reducible" stemming from the direct product of two intervals 
as one-dimensional Newton diagrams. In turn, the irreducible convex Newton dia- 
grams in R 2 are likely exhausted by different triangles, connected to the algebras 
sl 3 (R) and gl 2 (R) X R r+1 , r — 2, 3, 4, . . . (see Figures 4.1, 4.3). Those algebras play 
a role of primary algebras. At least, we could not find any other convex Newton 
diagrams in R 2 (see conjecture 1). 

Before description of our present knowledge about convex Newton diagrams for 
the R k case, let us recall the regular representation of the Lie algebra sZjv(R) 
given on the flag manifold which acts on the smooth functions of N(N — l)/2 real 
variables Zi i+ q , i = 1, 2, . . . , (N — 1), q — 1,2, ... , (N — i). The explicit formulas 
for generators in the Chevalley basis are given by (see pH ) 





i-l 






(126) 



:!8 
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where we use a notation ej, /,-, ft,-, i = 1,2,..., (iV— 1) for generators of positive and 
negative roots, and the Cartan generators, respectively. The algebra si/v(R) is gen- 
erated by generators e*, fi,hi. If rij are non-negative integers, the finite-dimensional 
irreducible representation of sZjv(R) will occur in the form of inhomogeneous poly- 
nomials in variables Zi^+ q . The highest weight vector is characterized by the integer 
numbers rij, i = 1, 2, . . . , (iV — 1). If some integers rii become zero, the representa- 
tion (126) is degenerated and acts on a space of dimension lower than N(N — l)/2 
(see, e.g., (88) as an example of a degenerate representation of s^R)). 

In connection with the space R 3 , at least two irreducible convex Newton 
diagrams appear: a tetrahedron, related to a degenerate finite-dimensional 
representation of s^4(R) (the highest weight vector is characterized by two zeroing 
and one non-zeroing components), and certain geometrical bodies, corresponding to 
the regular representation of s^R) on the flag manifold (see (126) at N — 3) and 
a representation of sj»2(R), respectively. We do not know whether those exhaust 
all irreducible Newton diagrams or not. 

For the general case R k , firstly, we wish to indicate, that there also exist Newton 
polyhedra related to the almost degenerate representation of slk+i(R) (in formulas 
(126) all integers except the one characterizing the highest weight vector become 
zero). Also there exist some convex geometrical figures related to some finite- 
dimensional degenerate and, sometimes, regular representations of sZfc + i_i(R), 
i > 0.0 It is easy to show that the algebra so/c+i(R) acts on R k in a degen- 
erate representation 

d 9 

Hij — Xi ~q^~ x j i * 3 i h 3 — 1,2,... ,fc 

fc d d 

Gi = ^2 x i x iJ> — + 0- + x f)^r~ ~ nxi , i = 1,2,... ,k 

j=i,M7 ■' '' (127) 

where Xi,i = 1,2, ... ,k are real variables and n is a number. If n is a non-negative 
integer, the algebra (127) possesses the finite-dimensional invariant sub-space, [^] 
which has the form of a polyhedron in R k and coincides with a finite-dimensional 
invariant sub-space of the algebra sifc+i(R) in the almost degenerate representation, 
acting on R k . Making a gauge transformation, any symmetric bilinear combination 
in the generators (127) of the algebra sofc+i(R) with real coefficients can be reduced 
to the form of a Laplace-Beltrami operator with some metric tensor plus a scalar 
function giving rise to exactly-solvable operators of the second type (for the proof 
sec [[l7| and the discussion in the end of Section 4.1). 

The above procedure can be extended to the case of the Lie super-algebras and 
also the quantum algebras. For the former case, our experience is very limited. The 
only statement one can make is that for linear operators acting on functions of one 
real and one Grassmann variable (or, equivalently, on two-component spinors in on 
real variable) there can appear finite-dimensional invariant sub-spaces other than 
finite-dimensional representations of the algebra osp(2,2) (see Section 3.1 and the 



15 In order to decide what algebras act on R k , it is necessary to perform some analysis of the 
Young tableau. I am grateful to W. Riihl for an enlighting discussion of this question. 

16 The corresponding finite-dimensional representation is unitary and reducible (see, e.g., the 
discussion in the end of Section 4.1 and also in Jl7t). 
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discussion in PI). Nevertheless, those sub-spaces are connected with the represen- 
tations of polynomial elements of the universal enveloping algebra of the algebra 
osp(2, 2) and, correspondingly, they can be constructed through finite-dimensional 
representations of the original super-algebra osp(2, 2). Probably there are no more 
irreducible convex Newton diagrams in this space. For the latter (quantum alge- 
bras), there exists an irreducible Newton diagram in R - an interval (see Section 
2) - but we could not find any irreducible (non-trivial) Newton diagrams in R k 
for k > 1; however, at least, in R 2 they occur, if a quantum plane is introduced, 
X V — 1V x i instead of an ordinary one, xy = yx, or a quantum hyperplane p0[ | for 
the case of more than two variables P]. 

It is worth noting that one can construct quasi-exactly-solvable and exactly- 
solvable operators considering "mixed" algebras: a direct sum of the Lie algebras, 
Lie superalgebras and the quantum algebras, realized in the first-order differential 
and/or difference operators, possessing finite-dimensional invariant subspaces. In 
this case the Newton diagrams are reducible. 

As for the cases considered in Sections 1-4, the problem of reducibility of 
general quasi-exactly-solvable and exactly-solvable operators to the Schroedinger- 
type operators is open. 

Above we have described the linear operators which possess a finite-dimensional 
invariant sub-space with a polynomial basis. Certainly, making changes of 
variables and gauge transformations one can obtain the operators with an invariant 
sub-space with a non-polynomial basis (but emerging from a polynomial one, see 
e.g. (11)). Those 'induced' linear spaces we consider as the equivalent to the orig- 
inal, polynomial ones. A general problem can be posed: is it possible to describe 
all linear differential operators possessing a finite-dimensional invariant sub-space 
with a basis in a certain explicit formp|, which cannot be reduced to a basis of poly- 
nomials using a change of variables and a gauge transformation? Equivalcntly, this 
problem can be re-formulated as a problem of a description of finite-dimensional lin- 
ear spaces of functions, where some linear differential operators can act. Since the 
original aim of the present investigation was mainly the construction of the quasi- 
exactly-solvable Schroedinger operators |l4|, [| [l6|, [l^], we did not touch on 
the general problem (except for a short remark on page 4) . A solution of this prob- 
lem plays a crucial role for a complete classification of the Schroedinger operators 
possessing a finite-dimensional invariant sub-space with an explicit basis in func- 
tions. Although it is rather surprising, it is very likely that this general and abstract 
problem has a constructive solution, at least for the one-dimensional Schroedinger 
operators. This work is in progress. 
Concluding remarks. 

(I). A crucial approach to the problem of classification of linear operators pos- 
sessing a finite-dimensional invariant sub-space(s) with a basis in functions in an 
explicit form|^] is connected to a theory of Riemann surfaces. 

17 I am grateful to O.V. Ogievetsky for an interesting discussion of this subject and valuable 
comment. 

18 It is quite difficult to define precisely the meaning of the words "space with a basis in a 
certain explicit form." One of the examples is the linear space < 1, x, /l(cc), fo(x), . . . , f n (x) >, 
where fi,i = 1,2, ...n are any functions on R 1 . Other examples of such spaces are the spaces (1), 
(11), (29), (77), (89), (113). 

19 Recently, a classification of linear differential operators possessing a finite-dimensional in- 
variant sub-space with a basis in monomials was completed P3], 

20 see footnote 18 
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Take a symmetric operator depending on a real parameter polynomially. and 
possessing the infinite discrete set of eigenvalues. Searching of eigenvalues is equiv- 
alent to solving the transcendental secular equation. It is rather natural that the 
roots of this equation (eigenvalues) form an infinite-sheet Riemann surface as a pa- 
rameter manifold. Generically, this surface is characterized by square-root branch 
points (see the excellent paper by Bender and Wu |^5| and also [^6[ p8|). Exis- 
tance of iV-dimensional invariant sub-space with an explicit basis implies, that the 
secular equation is reduced to two equations: an algebraic one and a transcenden- 
tal one. It reflects the fact, that the original determinant factors into the product 
of two determinants. Equivalently, the infinite-sheet Riemann surface is split off 
into two disjoint Riemann surfaces: A-sheet one and an infinite-sheet one. If there 
exist k finite-dimensional invariant sub-spaces (for instance, which are embedded 
consequently one into the other, or, another case, any two have no intersection), 
the original, infinite-sheet Riemann surface is split off into k finite-sheet Riemann 
surfaces and an infinite-sheet one. Namely, such a situation occurs for a general 
quasi-exactly-solvable operator. As for exactly-solvable operators, the correspond- 
ing infinite-sheet Riemann surface is split off completely for separate sheet sp^j Fol- 
lowing this idea, in [M |2^] a complete classification was done of second-order dif- 
ferential operators in R related with sZ2(R)-hidden algebra. This insight allowed 
to separate the quasi-exactly-solvable and exactly-solvable operators from generic 
operators, possessing abstract finite-dimansional invariant sub-spaces. I hope that 
such an approach can lead to an abstract classification of linear operators. 

(II) . The above-discussed questions about the existence of convex Newton dia- 
grams (equivalently, the finite-dimensional linear spaces with polynomial basis) as 
an invariant sub-space of some linear differential operator are related to a classifica- 
tion of finite-dimensional Lie algebras of first-order differential operators possessing 
an invariant sub-space. Originally, this problem was formulated by Sophus Lie, who 
solved this problem for the algebras on C 1 and gave a complete classification of the 
algebras of vector fields on C 2 (see |2^| and references therein). A complete classifi- 
cation of Lie algebras of first-order differential operators on C 2 has been presented 
just recently |l9|], while the same problem for the algebras on R 2 has not been 
solved so far, although the vector fields on R 2 have been described |22). The same 
problem for algebras on C 3 has also been discussed by S. Lie, who presented 35 al- 
gebras of vector fields, emphasizing that this list is not complete. As far as we know 
at present, similar questions about a classification of the Lie super-algebras, both 
vector fields and first-order differential operators in real and Grassmann variables, 
are never discussed in literature. The same situation takes place for the algebras 
of first-order finite-difference operators on the real line, or a quantum hyperplane 
HI- 
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APPENDIX 



In connection to the finite-dimensional representations of the algebras of first- 
order differential (difference) operators some operator identities occur [^| . 
1. The following operator identity holds: 



(J+)" +1 ^(x 2 a ;c -nxr+ 1 =.x 2 "+ 2 ^+ 1 , d x = ± n = 0,1,2, 



(A.l) 



Proof. The proof is straightforward: 

(i) the operator annihilates the space of all polynomials of power not 
higher than n, V n (x) = Spanlx 1 : < i < n}; 

(ii) in general, an (n + l)-th order linear differential operator annihilating V n {x) 
must have the form B(x)d" +1 , where B(x) is an arbitrary function and 

(iii) since is a graded operator, dcg(J+) = +1, p| deg((J+)" + ) =n+l, 
hence B(x) — bx 2n+2 , while clearly the constant 6 = 1. 

□ 

It is worth noting that taking the power in (A.l) different from (n+ 1), the l.h.s. in 
(A.l), generally, will contain all derivative terms from zero up to (n + l)-th order 
with non-vanishing coefficients. 

The identity (A.l) has a Lie-algebraic interpretation. The operator (</+) is the 
positive-root generator of the algebra of first-order differential operators (the 
other s/2-generators are J° = xd x — n/2 , J~ — d x , see (2)). Correspondingly, the 
space V n {x) is nothing but the (n + l)-dimcnsional irreducible representation of 
sl2- The identity (A.l) is a consequence of the fact that (Jn) n+1 = in matrix 
representation 

Another Lie-algebraic interpretation of (A.l) is connected with occurence of some 
relations between elements of the universal enveloping algebra of three-dimensional 
Heisenberg algebra, Hi: {P,Q, 1}- Let [P,Q] — 1, then 

(Q 2 P - nQ) n+1 = Q 2n+2 P n+1 , n = 0, 1, 2, . . . (A.2) 

and certainly this identity is more general than (A.l). Taking different represen- 
tations for P and Q in terms of differential operators (others than the standard 
one, P = d x and Q — x), one can get various families of operator identities, other 
than (A.l). Also the formula (A.2) has a meaning of a formula of an ordering: 
Q-operators are placed on the left, P-operators on the right. 

One can easily check, that once two operators obey [P, Q] — 1, then 

J+ = Q 2 P - nQ , 

J° = QP-n/2 , (A.3) 
J- = P, 

obey s^-algebra commutation relations. The representation (2) is a particular case 
of this representation. In fact, (A.3) is one of possible embeddings of the algebra 
s/2 into the universal enveloping algebra of three-dimensional Heisenberg algebra. 



22 so Jn maps x k to a multiple of x k+1 

23 It is a particular case of more general statement: in the algebra gl^ a positive (negative) 
-root generator J± a in a finite- dimensional representation of the dimension d, taken in power of 
of the dimension d is egual to zero, (J± a ) d = 0. 
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There exist other algebras of differential or finite-difference operators (in more 
than one variable), which admit a finite-dimensional representation. This leads to 
more general and remarkable operator identities. In fact, (A.l) is one representative 
of an infinite family of identities for differential and finite-difference operators. Also 
there occur generalizations of the relation (A. 2) for a certain polynomial elements 
of the universal enveloping algebra of (2n + l)-dimensional Heisenberg algebra. 

2. The Lie-algebraic interpretation presented above allows us to generalize (A.l) 
for the case of differential operators of several variables, taking appropriate powers 
of the highest-positive-root generators of (super) Lie algebras of first-order differen- 
tial operators, possessing a finite-dimensional invariant sub-space. First we consider 
the case of sl^ . There exists a representation of sl^ (R) as differential operators on 
R 2 (see Section 4.1, eq.(88)). One of the generators is 

•^K 71 ) = x 2 d x + xydy — nx . 

The space (89), V n {x,y) — Span{x i y : ' : < i + j < n} is a finite-dimensional 
representation for sl^, and due to the fact that (J\ (n)) n+1 = on the space V n {x, y), 
we have 

(4(n)) n+1 = (x 2 d x + xydy - nx) n+1 = 

^{ n+ k % 2n+2 - k y k d n x +1 - k d h y ■ (A.4) 

k=0 ^ ' 

This identity is valid for y € R (as described above), but also if y is a Grassmann 
variable, i.e. , y 2 — In the last case, J\{ri) is a generator T + of the algebra 
osp(2,2) (see (31)). 

More generally (using algebra slk+i instead of SI3), the following operator iden- 
tity holds: 



(j£ 1 {n)) n =(xi(^T x m d x 



n)) 



n+1 



c 



n+l 



r n+l \^ 

h+h+—+3k=n+l 



x l x 2 ■ ■ ■ x k u x 1 u x 2 ■ ■ ■ U x k 



(A.5) 



where C 1 ™^ . are the generalized binomial (multinomial) coefficients. If x £ 
R k , then J^ _1 (n) is a generator of the algebra sife+i(R), generated by D(fi) (see 
(126)). While some of the variables x's are Grassmann ones, the operator J^" 1 ^) 
is a generator of a certain super Lie algebra of first-order differential operators. 
The operator in the l.h.s. of (A.5) annihilates the linear space of polynomials 
V n {xi,X2, ■■■Xk) = Span{x 3 ^X2 . . . x J k k : < ji + j 2 + ■ ■ ■ + jk < n}, which is 
represented by tetrahedron Newton diagram in R k . 
Considering (2k + l)-dimensional Heisenberg algebra: 

[Pi,Qj] = Sij, [Pi, = 0, [Q u Qj] = 0, i,j = 1,2, . . .k 



In this case just two terms in the l.h.s. of (A.4) survive. 
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one can rewrite (A. 5) in more general form 

k 

(Ql(J2 QmPm - n)) n+1 = 
m—1 

Q n + 1 £ C n+l^ Q , Q n Q^phph ph j (A . 6) 

jx+h+— +Jfc=n+1 

(cf.(A.2)). It is clear, the formulas for sifc + i-generators at k > 1 (see (126)) in 
terms of Pi,Qi, analogous to (A. 3), can be easily derived. 

3. For the case of the algebra gfoCR) x R r+1 , r = 1,2,... (see (112)) the situation 
becomes rather sophisticated 

(J*) n+1 = (x 2 d x + rxyd y - nx) n+1 = 
n+l / _|_ 1\ 

V r k [ n I x 2n + 2 - k y n+1 - k d k ^+ 1 - k + Q rtn , (A.7) 
k=o V k ' 
where Q\ %n = and also 

at n=0 

r(r — l)x 2 yd y at n=l 

r(r - l)x 3 y[3ryd 2 + 3xd 2 y + (r ~ 2)d y ] at n=2 

(cf. (A.4)). For arbitrary n and r 

r(r — V)n[n + 1) 

2 

[s* 1 ^- 1 ^, + r(n - lJa^-V^ -2 ^ + - - 2)x 2 "-V#T 3 ^ + . . .] 

+ [Zoiuer order derivatives]. 

The r.h.s. of (A.7) still has quite regular structure, which can be analyzed using 
the following trick: variables j/'s in the l.h.s. of (A.4), (A.7) can be related to each 
other through the very simple connection 

V(A.l) = y\ A A) > 

and then J\ (see (A.4)) coincides with J*. 

4. The above-described family of operator identities can be generalized for the 
case of finite-difference operators with the Jackson symbol, D x (see, e.g., |ll|), 
obeyed the following Leibnitz rule 



f {x) - f(q 2 x ) 
(l-q 2 )x 



D x f(x) = J \ ; + f(q 2 x)D x , 



instead of the ordinary, continuous derivativej^j Here, q is an arbitrary complex 
number. The following operator identity holds: 

(J+)" +1 = (x 2 D x - {n}x) n+1 = q 2n ^x 2n+2 D^+\n = 0, 1, 2, . . . 

(A.8) 

i _ 2n 

(cf. (A.l)), where {n} — -j^r is a so-called g-number. The operator in the l.h.s. 
annihilates the space (1), V n (x). The proof is similar to the proof of identity (A.l). 



25 The definition of the Jackson symbol used here is slightly different than that presented ir 
Section 2. 
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From the algebraic point of view the operator J+ is the generator of a q-deformed 
algebra s^(R-)q of first-order finite-difference operators on the line: 
J® = xD — h, J~ = D, where n = ^|2k+2}^ ( see @ anc ^ a ^ so HI)' obeying 
the commutation relations (22). 

One can show, that once two operators obey PQ — q 2 QP = 1, then 

{Q 2 P - {n}Q) n+1 = g 2n(n+l)Q2n+2pn+l ; n = , 1, 2, . . . (A.9) 

holds and 

J + = Q 2 P - {n}Q , 

J° = QP-n , (A.10) 
J" = P, 

obey s^2(R)q-algebra commutation relations (22). It is evident, that this represen- 
tation is more general than the representation (21). 

An attempt to generalize (A. 4) replacing continuous derivatives by Jackson 
symbols immediately leads to requirement to introduce the quantum plane and 
g-differential calculus ||(| 

xy = qyx , 

D x x = 1 + q 2 xD x + (q 2 - l)yD y , D x y = qyD x , 
D y x = qxDy , D y y = 1 + q 2 yD v , 

D x D y = q~ 1 D y D x . (A.ll) 
The formulae analogous to (A. 4) have the form 

(J 2 1 (n))" +1 = (x 2 D x + xyD y - {n}x) n+1 = 

"+1 / , -i \ 

V- q 2n*-n(k-2) + k(k-l) (n + x 2n+2- ky k D n+l^k D k _ {A U) 

fc=0 \ k Jq 

where 

are q-binomial coefficient and q-factorial, respectively. As in all previous cases, if 
y E R, the operator J^in) is one of generators of the q-deformed algebra SI3 (R) q 
of finite-difference operators, acting on the quantum plane and having the linear 
space V n [x, y) — Span{x l y J : < i + j < n} as a finite-dimensional representation; 
the l.h.s. of (A. 12) annihilates T' n (x,y). If y is a Grassmann variable, J 2 (ri) is 
a generator of the q-deformed superalgebra osp(2, 2) q possessing finite-dimensional 
representation. 

Introducing two couples of the operators Pi,2,Qi,2 obeying the following rela- 
tions: 

Q1Q2 = 17Q2Q1 , 

P1Q1 = 1 + q 2 QiPi + (q 2 - i)Q 2 p 2 , P1Q2 = qQ 2 Pi , 

P2Q1 = qQlP2 , P2Q2 = 1 + q 2 Q2P2 , 

P 1 P 2 =q- 1 P 2 P 1 . (A.13) 
It is easy to derive a family of abstract identities similar to (A. 6) at k = 2 
{Q\Pi + Q1Q2P2 - {n}Qi) I1+1 = 
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n+1 / , i \ 

y- (7 2r 1 2 -„(fe-2)+/c(fc-l)/«+ i \ Q 2n+2-k Q kpn+l-kpk > (A U) 

fc=0 \ k Jq 

while (A. 12) is a particular case corresponding a special choice of a representation 
(A.13) in the form of (A.ll). 

Introducing a quantum hyperplane [ pp| , one can generalize the whole family 
of operator identities (A. 5) replacing first of all continuous derivatives by finite- 
difference operators and then introducing the abstract operators like it was done 
above for the case of (A.11)-(A.12). 

The Lie-algebraic interpretation of above operator identities allows to make a 
general conclusion that an existence of finite- dimensional representations of the Lie 
algebra gl^ of differential operators leads to an appearence of some specific operator 
identities analogous to those families described above. 
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